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The literature on self-diffusion of polymers in solution, and on tracer diffusion of probe polymers 
through solutions of matrix polymers, is systematically reviewed. Virtually the entirety of the 
published experimental data on the concentration dependence of polymer self — and probe-diffusion 
is represented well by the same functional form. This form is the stretched exponential exp(— ac"), 
where c is polymer concentration, a is a scaling prefactor, and v is a scaling exponent. 

Correlations of the scaling parameters with polymer molecular weight, concentration, and size are 
examined, and compared with predictions based on the form's hydrodynamic and renormalization- 
group derivations, a increases markedly with polymer molecular weight, namely a ~ M x for x « 1. 
v is ~ 0.5 for large polymers (M larger than 400 kDa or so), but increases toward 1.0 or so at 
smaller M. Scaling parameters for the diffusion of star polymers do not differ significantly from 
scaling parameters for the diffusion of linear chains of equal size. 
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I. INTRODUCTION 

The objective of this review is to present the experi- 
mental phenomenology for polymer self-diffusion and for 
diffusion of tracer polymers through polymer solutions. 
We proceed by examining how the diffusion coefficient 
depends on polymer concentration c, polymer molecular 
weight M, solvent quality, and other variables. We exam- 
ine how rapidly a probe polymer of molecular weight P 
diffuses through a matrix polymer of molecular weight 
M, perhaps as the concentrations of the matrix and 
probe polymers are varied. After establishing the phe- 
nomenological behavior, we compare the phenomenology 
with various classes of theoretical model. 

This is not a review of the extremely extensive theoret- 
ical literature on polymer diffusion in solution. For such 
reviews, note works of Graesslej*i, Tirrelfi, Pearson^, 
Skolnick 4 , Lodge 5 , and (more recently but less di- 
rectly) McLeish&. Recent papers by Schweizer and 
collaborators^^ include extensive references to the more 
recent literature. We will note classes of models, but not 
their underlying derivations. We also do not consider 
melt systems, or polymer mutual diffusion. 

In comparing phenomenology with model predictions, 
it will remain critical to distinguish between properties 
that are consistent with a particular model, but do not 
actually prove it, and properties that require or refute 
the correctness of a particular model. For example, 
SkolnichiS has shown that the nominal signature of rep- 
tation D s ~ M™ 2 is found in computer models of poly- 
mers, in a model system in which the chains are very cer- 
tainly not reptating. A finding D s ~ M~ 2 is thus consis- 
tent with many tube-type solution models, but does not 
prove their correctness. On the other hand, most tube- 
type models require that polymer chains move via con- 
strained diffusion for shorter times. In these models, the 
shorter-time mean-square displacement of a single chain 
must scale more slowly than linearly with time. If the 
shorter-time mean-square displacement were instead lin- 
ear in elapsed time, many constrained-diffusion models 
would be rejected by experiment. 



Section II presents a theoretical background for the 
experimental papers reviewed here. Sections III and IV 
treat, respectively, data on (i) self-diffusion coefficients of 
polymers of solution and (ii) data on probe diffusion co- 
efficients in polymer matrix solutions. Section V briefly 
remarks upon other papers on polymer self- and probe- 
diffusion that do not lend themselves to the analytic 
approach applied here. A systematic study of the phe- 
nomenological parameters obtained by our analysis ap- 
pears in Section VI. Section VII summarizes conclusions. 
Tables of fitting parameters appear as Appendices. 



II. THEORETICAL BACKGROUND 

This Section presents a consistent nomenclature for 
polymer diffusion coefficients and describes methods 
whereby which such diffusion coefficients are measured. 
A short description is given of the literature on diffu- 
sion in multimacrocomponcnt solutions. The literature in 
question provides the fundamental basis for interpreting 
the experimental measurements. Finally, phenomenolog- 
ical classes of models for polymer dynamics are identified. 



A. Nomenclature for Diffusion Coefficients 

This section sets out a consistent nomenclature for the 
diffusion coefficients that we are reviewing. In a solution 
containing a solvent and one macromolecular species, two 
physically distinct diffusion coefficients usefully charac- 
terize macromolecule diffusion. One of these, the two- 
particle or mutual diffusion coefficient D m , describes via 
Fick's law 



J = D m Vc 



(1) 



the relaxation of a concentration gradient. Here J is the 
diffusion current and c is the local instantaneous macro- 
molecule concentration. The other diffusion coefficient, 
the single-particle or self diffusion coefficient D s , de- 
scribes the diffusion of a single macromolecule through a 
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uniform solution of elsewise identical macromolecules. In 
a simple macromolecule solution, quasi-elastic light scat- 
tering spectroscopy (QELSS) measures the mutual diffu- 
sion coefficient of the macromoleculesii*i£. Pulsed-field- 
gradient nuclear magnetic resonance (PFGNMR) mea- 
sures the self-diffusion coefficient. 

Solutions containing a solvent and two macromolecule 
species show a more complicated diffusive behavior. The- 
oretical treatments of this issue are described in the next 
section. In general, in two-macrocomponent solutions 
the temporal evolution of transient concentration fluc- 
tuations is described by two relaxation times, each de- 
scribing the relaxation of a coupled mode involving the 
concentrations of both macromolecule species. An in- 
teresting special case arises if one macrocomponent, the 
probe, is adequately dilute, while the other macrocompo- 
nent, the matrix, may be either dilute or concentrated. In 
this special case, the diffusion of the probe species is gov- 
erned by a single-particle diffusion coefficient, namely the 
probe diffusion coefficient D p . In the special case that the 
matrix species and solvent are isorefractive, so that the 
index of refraction is independent of the relative amounts 
of matrix polymer and solvent, the matrix species scat- 
ters next to no light. If in this special case light scattering 
by a dilute probe species dominates scattering by the so- 
lution (which further requires that solvent scattering be 
sufficiently weak), the diffusion coefficient measured by 
QELSS ia 12 i 13 the probe diffusion coefficient D p . 

A variety of physical techniques have been used to mea- 
sure probe diffusion through a polymer solution. In the 
fluorescence recovery after photobleaching (FRAP) tech- 
nique, a small fluorescent label is attached to the probes. 
An intense pulse of light (the 'pump' laser beam) is then 
used to destroy ( "bleach" ) the fluorescent labels in some 
regions of the solution, the regions sometimes being de- 
fined via a holographic grating. Much weaker laser il- 
lumination (the 'probe' beam) is then used to monitor 
the recovery of the fluorescence intensity as unbleached 
fluorescent molecules diffuse back into the regions (s) in 
which the bleaching occurred. 

An alternative to FRAP is Forced Rayleigh Scatter- 
ing (FRS). In an FRS experiment, an intense laser pulse 
is used to generate a holographic brightness grating in 
solution. The grating selectively alters a photosensitive 
part of the probe molecule, thereby creating an index of 
refraction grating in the solution. A much weaker probe 
beam is then used to monitor the diffusive relaxation of 
the induced grating. Significant complications may arise 
if the photomodified and non-photomodified forms of the 
probe species differ significantly in their diffusive proper- 
ties. If the probe species is dilute, with either FRAP or 
FRS the time dependence of the recovery profile is deter- 
mined by the single-particle (probe-) diffusion coefficient 
of the probe molecules^. 

A physically-distinct alternative to FRAP and FRS 
is fluorescence correlation spectroscopy (FCS), in which 
optical methods are used to measure fluctuations in the 
number of fluorescent molecules in a small volume of 



space. In nondilute solution the diffusion coefficient 
measured by fluorescence correlation spectroscopy (FCS) 
varies from D p to D m as the fraction of macromolecules 
that bear fluorescent tags is varied from small to largeii 

When is a probe species dilute? In general, the care- 
ful experimenter makes an adequate control study of the 
effect of varying the probe concentration c p . D p is the 
extrapolation of the measured probe diffusion coefficient 
to zero probe concentration. In some systems D p is 
substantially independent of c p , and no extrapolation is 
needed. In other systems (see below) D p depends sig- 
nificantly on c p , and extrapolation to c p — > must be 
performed. Several authors have measured the initial 
slope kd of the dependence of D p on probe concentra- 
tion, and the effect of matrix concentration on kd- It 
is also possible to study-most published work has used 
QELSS-the relaxation of concentration fluctuations in a 
ternary solvent: macromolecule: macromolecule system in 
which neither macrocomponent is dilute. As reviewed 
below, such studies give information on the diagonal — 
and cross-diffusion coefficients of the two species and on 
their thermodynamic interactions. 

Strictly speaking, if FRAP or FRS is used to study 
the diffusion of a dilute labelled polymer through a 
matrix polymer solution, the labelled polymer differing 
from the matrix in that a fluorescent small-molecule la- 
bel has been attached, the measured diffusion coefficient 
should be termed a probe diffusion coefficient, because 
the molecules being tracked are chemically distinguish- 
able from the molecules not being tracked. However, for 
the purposes of this review, we distinguish between (i) 
the case in which the probe and matrix polymers are sub- 
stantially distinct in molecular weight or chemical nature, 
and (ii) the case in which the probe and matrix polymers 
are very nearly the same except for the presence of the 
label. The phrase 'probe diffusion coefficient' is reserved 
for case (i), while measurements under case (ii) are here 
treated as determinations of the self diffusion coefficient. 

The term tracer diffusion coefficient is synonymous 
with single-particle diffusion coefficient, and includes 
both the self and probe diffusion coefficients as special 
cases. The interdiffusion and cooperative diffusion co- 
efficients characterize the relaxation times in a ternary 
system in which neither macrocomponent is dilute. In 
the following the latter two terms are almost never used. 



B. General Theory of Diffusion in 
Multimacrocomponent Solutions 

What is the basis for using spectroscopic methods such 
as QELSS, FRAP, FRS, or FCS to measure D s or D p l 
In each case, there is a physical theory that links the cor- 
relation function determined spectroscopically to fluctu- 
ations in microscopic variables that describe the liquid. 
For example, for FCS the fluctuations in the fluorescent 
light intensity If (t) arise from fluctuations in the number 
of fluorescent molecules in the scattering volume. The 
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time correlation function C/(r) = (If(t)If(t +r)) of the 
fluorescent intensity follows faithfully the time correla- 
tion function C n {t) = (N(t)N(t + r)} of the (fluctuat- 
ing) number N(t) of labelled molecules in the scattering 
volume, so to determine the time dependence of C/(r) 
one only needs to determine Cjv(i) and vice versa. 

For QELSS, the spectrum is determined^ by the field 
correlation function (q, t) of the scattered light, which 
is in turn determined by the motion of the scattering 
molecules via 

g {1) (q,T) = / ]T ctiOtjexp[iq- (r;(t + r) -ry (<:))] \ , 

(2) 

where i and j label two of the N scatterers, Yi{t) is the 
position of scatterer i at time t, aj is the scattering cross- 
section of particle j, and q is the scattering vector. For 
the special case in which the scattering particles are di- 
lute (non-scattering particles do not need to be dilute), 
correlations between the positions of two distinct scat- 
terers at two times vanish. The correlation function for 
probe scattering then reduces to 

^ 1) fer) = /^ a ?exp[iq.Ar i (r)]^, (3) 

in which Ar^(r) = Ti(t + r) — Yi(t)). 

There are three general approaches to calculating 
g P l \q, r). First, scatterers can be treated as objects hav- 
ing hydrodynamic and direct, e.g., excluded volume, in- 
teractions, and the effect of the intermacromolecular in- 
teractions on dgp^ (q, r) / dr can be calculated. Second, it 

can be recognized that J^f ex P[*q - ( r i(*)] is the q th spatial 
Fourier component of the scatterer concentration at time 
t, and semicontinuum hydrodynamic methods and the 
Onsager regression hypothesis can be used calculate the 
average temporal evolution of a fluctuation. Third, it can 
be proposed that the scatterers perform simple Brown- 
ian motion, with the probability distribution function for 
Ar^(r) being applied to calculate gp 1 ' (q, t). 

Light scattering spectra of non-ideal single-component 
macromolecule solutions, including direct interac- 
tions and hydrodynamic interactions at the Oseen 
level, were initially calculated by Altcnberger and 
Deutsche. Calculations of light scattering spectra 
of non-dilute many-component macromolecule solutions 
soon followedii*i2ii&. In the latter papers, diffusion coeffi- 
cients were modeled on the lines of Kirkwood, et al— , the 
diffusion coefficients being written as products of thermo- 
dynamic derivatives (dpn/dcj) obtained from the inter- 
macromolecular forces, and phenomenological transport 
coefficients fiy (here i and j label chemical species). If 
neither macromolecular species was dilute, the light scat- 
tering spectrum obtained by QELSS is predicte d 11 ' 16 to 
contain two relaxation modes. Even if only one of the 
two macromolecule species scattered light, both modes 
are predicted to be visible, in general, in the light scat- 
tering spectrum— 



A subsequent paper— extended these findings to the 
special case explored e.g., by fluorescence correlation 
spectroscopy, in which some macromolecules are tagged 
so that their movements can be tracked. For the 
case that the tagged macromolecules are dilute, it was 
shown— that fluorescence correlation spectroscopy and 
other equivalent techniques measure the self-diffusion co- 
efficient of the tagged probe macromolecules diffusing 
through the unlabelled matrix macromolecules. 

A Smoluchowski approach was used by Jonesi^ to ex- 
amine interacting spherical polymers, including the spe- 
cial case that one of the polymer species is dilute and 
tagged so that its motions could be observed. In the spe- 
cial case, the measured diffusion coefficient was predicted 
by Jones to be determined only by the hydrodynamic 
interactions of the tracer polymers with their (predom- 
inantly matrix-polymer) neighbors. While Jones does 
not use the same nomenclature, the diffusion coefficient 
of the tracers calculated in ref. may be recognized 
as the single-particle diffusion coefficient. The hydrody- 
namic approach culminates in analyses of Carter, et al- 
and Phillies 2 ^ of mutual and tracer diffusion coefficients, 
including full hydrodynamic and direct interactions be- 
tween interacting diffusing particles and reference frame 
issues. 

Extensive studies of diffusion of interacting polymers 
using semicontinuum and related means were made by 
Akcasu, Benmouna, Cohen and others. In 1987, Ben- 
mouna, et al^ 1 calculated dynamic scattering from a 
solution containing two polymer species, including ex- 
cluded volume terms with a Flory interaction parameter 
but ignoring hydrodynamic interactions between poly- 
mer chains. Hydrodynamic interactions were neglected 
in the sense that cross-species transport coefficients were 
assumed to vanish. The assumption that the hydrody- 
namic parts fly of the cross-diffusion tensor vanish can- 
not be true simultaneously in the solvent- and volume- 
fixed reference frames^. Polymer solutions can be very 
concentrated, so neglect of reference frame issues in en- 
tirely correct calculations may lead to complications in 
the physical interpretation of an assumption that hydro- 
dynamics has been neglected. 

Benmouna, et al. 2 — examined the special case of 
two species identical except for their optical scatter- 
ing cross-sections, showing results consistent with those 
of PhilliesiSiH, notably: If the matrix species scatters 
weakly, an interdiffusion mode describing single-chain 
motion is dominant at low probe concentrations. There is 
a thermodynamic regime relative to the spinodal where 
the diffusion equations describe phase separation. The 
Benmouna model was then extended to treat solutions of 
copolymers 2 ^, still with neglect of hydrodynamic interac- 
tions, and to treat homopolymer: copolymer mixtures 23 . 

In the same period, Foley and CohenSi analyzed con- 
centration fluctuations in polymer: polymer: solvent 
mixtures, using an ornate Flory-Huggins form for the 
thermodynamic free energy of the mixture, again neglect- 
ing interchain hydrodynamic interactions. The case of a 
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solvent isorefractive with the matrix polymer, in the pres- 
ence of a dilute scattering species, was treated. Foley 
and Cohen also examined systems in which both poly- 
mer species were nondilute, predicting that in this case 
the relaxation spectrum is characterized by three distinct 
relaxation times. In contrast, other calculations on sim- 
ilar models predict only two relaxation times. 

Roby and JoannyS^ improved the model of Benmouna, 
et al^i by incorporating interchain hydrodynamic inter- 
actions and by improving the model for direct chain- 
chain interactions. At elevated concentrations, the cor- 
rectness of the reptation model was assumed. The effect 
of reptation dynamics on concentration fluctuations was 
estimated with the approximation that a solution is ef- 
fectively a polymer melt in which mesoscopic polymer- 
solvent blobs play the role taken in polymer melts by 
monomer units. For simplicity, the model was restricted 
to systems containing equal amounts of two species hav- 
ing the same molecular weight. These restrictions are 
nontrivial to remove but exclude analysis of tagged-tracer 
experiments. A model calculation incorporating a simi- 
lar picture, but without the restrictions, was reported by 
HammoudaSii. 

An alternative analysis of a ternary solution of two 
polymers and a solvent was presented by Wang2£, who 
followed much of the earlier work on polymer solutions in 
assuming that the cross terms in the mobility matrix van- 
ish. Wang systematically analysed a variety of general 
and special cases, showing a range of interesting parame- 
ters that can be determined from light scattering spectra 
if the accuracy of his model is assumed. In addition to 
the tracer case (one polymer:solvent pair isorefractive, 
visible polymer dilute), Wang analysed the special case 
"zero average contrast" , in which the polymer refractive 
increments are of opposite sign, so that concentration 
fluctuations, that change the total polymer concentration 
without changing the local polymer composition, scatter 
no light. Wang showed that the QELSS spectrum of a 
zero-average-contrast system is almost always bimodal, 
though it may happen that one of the modes is much 
weaker than the other. 

Treatments of the field correlation function based on 
a Brownian motion description can be traced back to 
Berne and Pecora— , who treat the light scattering spec- 
trum of a solution of dilute, noninteracting Brownian 
particles. The motion of such particles is described as a 
series of random, uncorrelated steps, in which case from 
the Central Limit Theorem the probability distribution 
for particle displacements is 

G S (AR,<) = 



f <(AR) 2 > 



-3/2 



exp [-3(Ai?) 2 /2((AR) 2 >] , (4) 



where the mean-square particle displacement is related 
to the diffusion coefficient by 

<(AR) 2 ) = 6Dt. (5) 



Combining eqsEJand^ 

gp\q, t) ~ J dAr exp[zq • Ar] 



2tt 



((AR) 5 



-3/2 



exp[-3(Ai?) 2 /2((AR) 2 }], (6) 



which using eq [5] leads to 



(7) 



For Brownian particles in low-viscosity small-molecule 
solvents, the Stokes-Einstein equation 



D 



k B T 
6t:i]R ' 



(8) 



in which k B is Boltzmann's constant, T is the absolute 
temperature, n is the solution viscosity, and R is the 
particle radius, generally predicts accurately the particle 
diffusion coefficient. For diffusion in polymer solutions, 
two ad hoc extensions of this form are encountered: 

First, it might be the case that D is not predicted 
accurately by eq|H| In this case, one could formally define 
a microviscosity rj^ as 



6irDR' 



(9) 



and compares 77^ with the macroscopically-measured vis- 
cosities 77 of the solution and r\ Q of the solvent. The micro- 
viscosity is more commonly applied to describe diffusion 
of mesoscopic globular probe particles 2 ^, rather than to 
treat the diffusion of random-coil polymers in solution. 

Second, light scattering and other relaxation spectra 
are not always single exponentials. No matter what func- 
tional form gp\q,r) has, one may formally write 



ff^(«,T)~exp(-«r , D(r)T) 



(10) 



as the definition of D(r). In this equation, D(t) formally 
appears to be a time-dependent diffusion coefficient, 
which equally formally defines a frequency-dependent dif- 
fusion coefficient such as 

/>oc 

D{uj) = dr &qp(-itOT)D(T), (11) 
Jo 

and via several slightly different paths a frequency- 
dependent microviscosity 



k B T 



QnD{uo)R' 



(12) 



The real and imaginary parts of rj^{uj) can be brought 
into correspondence with the storage and loss moduli 
G'(lu) and G"{u). 

The second extension has serious physical difficulties. 
In particular, comparing eqs[3EI and llUI and noting from 
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reflection symmetry that terms odd in q vanish, one finds 
that the extension assumes that 

(exph^Ar) 2 ]) = exp[-g 2 <(Ar) 2 )]. (13) 

Equation 1131 would be correct if sequential random 
changes in Ylj=i exp[iq • rj] were described by a Gaus- 
sian random process, because in that case (and only in 
that case) the average on the lhs of eq EH would be en- 
tirely described by the mean-square average displacement 
via the expression ^(Ar) 2 ) seen on the rhs of this equa- 
tion. Brownian motion, in which particle displacements 
are described by the Langevin equation, is an example 
of a dynamics that generates a Gaussian random process 
for which eq^Dis correct. For extensive details, see ref 


Omitted from, but critical to, the above discussion 
are the consequences of Doob's First Theorem 2 ^. Doob 
treated the joint expectation value-what we would now 
call the correlation function-of random variables includ- 
ing those following Langevin dynamics. If ea 1131 is cor- 
rect, it is an inescapable consequence of Doob's Theo- 
rem that the relaxation spectrum is a single exponential. 
Conversely, if the spectrum is not a single exponential, 
then from Doob's theorem 

3 ( 1 )(q,r)^exp[- (Z 2 ((Ar) 2 }]. (14) 

Berne and Pecora's texliM is sometimes incorrectly 
cited as asserting that eq^|is uniformly correct for light- 
scattering spectra. The analysis in Berne and Pecora 1 -, 
which correctly obtains ea !13l refers only to a special-case 
system. In the system that these authors correctly ana- 
lyzed, particle displacements are governed by the simple 
Langevin equation, and particle displacements in succes- 
sive moments are uncorrelated. That is, this analysis 
refers to systems in which particle motion is observed 
only over times much longer than any viscoelastic relax- 
ation times. 

If particle motion were observed over times shorter 
than the viscoelastic relaxation times, which in a con- 
centrated polymer solution might be 1 second or more: 
Particle displacements in successive moments would 
be correlated. Equation ^5] would not be correct. 
logf^W (q, r)]/q 2 would not be proportional to the mean- 
square particle displacement during r. The path from eq 
12 to ea 1121 might be an interesting heuristic, but would 
not be consistent with the properties of particles execut- 
ing Brownian motion. Of course, with respect to polymer 
dynamics the interest in eas llOMTSI has been exactly the 
study of particle motions at short times, during which 
viscoelastic effects are apparent in D(t), but these are 
precisely the circumstances under which eq^|is incor- 
rect. 

In addition to the above, there are analyses of melt 
systems such as the very interesting work of Akcasu and 
collaborators^}^. Melt systems are not considered here. 

C. Phenomenological Forms for 



Comparison with Experiment 

The approach here is to compare experimental mea- 
surements of D s and D p with the functional forms and c, 
P, and M dependences predicted by various treatments 
of polymer dynamics. There are a very large number 
of proposed models. Most models fall into two major 
phenomenological classes, distinguished by the functional 
forms taken to describe D s (c). This section sketches pre- 
dictions of classes of model, not including their underly- 
ing physical rationales, in preparation for the compari- 
son. 

(1) In scaling-law models 32 , the relationship between 
D s , D v , and polymer properties is described by scaling 
laws such as 

D s = D x M^c- x , (15) 

where here 7 and x are scaling exponents, and D\ is 
a scaling prefactor, namely the nominal diffusion coef- 
ficient at unit molecular weight and concentration. In 
some cases, scaling laws are proposed to be true only over 
some range of their variables, or only to be true asymp- 
totically in some limit. On moving away from the limit, 
corrections to scaling then arise. For melts, some models 
derive a scaling law for D S (M) from model dynamics, and 
then predict numerical values for 7. (In melts, D s has 
no concentration dependence because c is constant.) For 
polymer solutions, more typically a scaling-law form is 
postulated. The theoretical objective is then to calculate 
the exponents 7 and x. 

Many scaling-type models propose a transition in so- 
lution behavior between a lower-concentration dilute 
regime and a higher-concentration semidilute regime. 
Scaling arguments do not usually supply numerical co- 
efficients, so there is no guarantee that an interesting 
transition actually occurs at unit value of a predicted 
transition concentration ct rather than at, say, 2ct. Cor- 
respondingly, the observation that a transition is found 
at 2ct rather than ct is generally in no sense a disproof 
of a scaling model, because in most cases scaling models 
do not supply numerical prefactors adequate to make a 
disproof. (Some level of rationality must be preserved. If 
a physical model leads to c t as the transition concentra- 
tion, and the nominally corresponding transition is only 
observed for 30-150 Ct, and then only in some systems, 
one is entitled to question if the observed transition cor- 
responds to the transition implied by the model.) 

Two transition concentrations are often identified in 
the literature. The first transition concentration is the 
overlap concentration c* , formally defined as the concen- 
tration c* = N/V at which 4nR 3 g N/(3V) = 1. Here N is 
the number of macromolecules in a solution having vol- 
ume V and R g is the macromolecule radius of gyration. 
In many cases, c* is obtained from the intrinsic viscos- 
ity via c* — n/[rj\ for some n in the range 1-4. The 
second transition concentration is the entanglement con- 
centration c e . In some papers, the entanglement con- 
centration is obtained from a log-log plot of viscosity 
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against concentration by extrapolating an assumed low- 
concentration linear behavior and an assumed higher- 
concentration power-law behavior (c x for, e.g., x = 4) 
to an intermediate concentration at which the two forms 
predict the same viscosity; this intermediate concentra- 
tion is taken to be c e . In other papers, the entangle- 
ment concentration is inferred from the behavior of the 
viscoelastic moduli; in particular, an onset of viscous re- 
covery in the melt or solution is taken to mark c e . 

(2) In exponential models^^i, the concentration de- 
pendence of D s becomes an exponential or stretched ex- 
ponential in concentration 

D s = D exp(-ac u ). (16) 

Here D a is the diffusion coefficient in the limit of infinite 
dilution of the polymer, a is a scaling prefactor, and v is 
a scaling exponent; v — 1 for simple exponentials. Un- 
der the circumstance that the probe and matrix molec- 
ular weights P and M differ, an elaborated form of the 
stretched exponential is 

D p = D Q p- a exp(-ac"P~<M s ), (17) 

where 7, a, and S are additional scaling exponents, D Q 
now represents the diffusion coefficient in the limit of 
zero matrix concentration of a hypothetical probe poly- 
mer having unit molecular weight, and P~ a describes 
the dependence, on probe molecular weight, of the dif- 
fusion coefficient of a dilute probe molecule. On setting 
a = 7 = and freeing D Ql ea. 1171 becomes a parameteri- 
zation of the matrix molecular weight dependence of D p 
for a fixed probe molecular weight. 

In the derivationa 33 i 34 i 35 i 36 of the stretched- 
cxponential models, functional forms and numerical 
values for exponents and pre-factors were both obtained, 
subject to various approximations. The latter two 
derivations assume that chain motion is adequately 
approximated by whole-body translation and rotation, 
which may be appropriate if P m M, but which is not 
obviously appropriate if P and M are substantially 
unequal, since in this case whole-body motion of one 
chain and local modes of the other chain occur over the 
same distance scale. 

Some exponential models^LiS also include a transi- 
tion concentration, namely a transition between a lower- 
concentration regime in which some transport coefficients 
show stretched-exponential concentration dependences 
and a higher-concentration regime in which the same 
transport coefficients show power-law concentration de- 
pendences. The transition concentration, which experi- 
mentally is sharply defined2LiS, is here denoted c + . The 
lower-concentration regime is the solutionlike regime; the 
higher-concentration regime is the meltlike regime. Such 
transitions are seen in some but not all viscosity data—, 
generally but not always 39 at very high concentrations 
c[rj\ > 35. A solutionlikc-mcltlikc transition is very rarely 
apparent in measurements of D s or D p . 

There has been interest in derivations of concentra- 
tion and other dependences of transport coefficients from 



renormalization group approaches. Power-law and expo- 
nential forms caniS both follow from a renormalization- 
group approach, depending on the location of the sup- 
porting fixed point. The stretched-exponential form 
is££ an invariant of the Altenberger-Dahler— Positive- 
Function Renormalization Group. 

Our analysis will examine whether either of these func- 
tions describe experiment. While a power law and a 
stretched exponential both can represent a range of mea- 
surements, on a log-log plot a power law is always seen 
as a straight line, while a stretched exponential is always 
seen as smooth curve of nonzero curvature. Neither form 
can fit well data that is described well by the other form, 
except in the sense that in real measurements with ex- 
perimental scatter a data set that is described well by 
either function is tangentially approximated over a nar- 
row region by the other function. 

It would also be possible to divide theories of poly- 
mer dynamics into classes based on assumptions as to 
the nature of the dominant forces in solution. Asser- 
tions as to the dominant forces are independent of asser- 
tions as to the functional form used to describe D s (c). 
The major forces common to all polymer solutions are 
the excluded-volume force that prevents polymer chains 
from interpenetrating each other, and the hydrodynamic 
force that creates correlations in the displacements of 
nearby chains. In some models, excluded-volume forces 
(topological constraints) are assumed to dominate, hy- 
drodynamic interactions serving primarily to dress the 
monomer diffusion coefficient. In other models, hydro- 
dynamic interactions between nearly chains are assumed 
to dominate, while chain-crossing constraints provide at 
most secondary corrections. 

In addition to the generic forces, chains have 
chemically-specific interactions including in various cases 
van der Waals, hydrogen bonding, and electrostatic 
forces. These interactions substantially modulate the 
properties of particular polymers. However, diffusion 
coefficients and viscoelastic parameters of most neutral 
polymers show highly-characteristic "polymeric" behav- 
iors, almost independent of the chemical identity of the 
polymer, implying that general polymer properties do 
not arise from chemically-specific interactions. 

In the following, eas. H5tfr7l will be systematically com- 
pared with the literature on D s and D p . The follow- 
ing largely treats diffusion by neutral polymers in good 
and theta solvents. There is rather little data on self- 
diffusion of random coil polyelectrolytes. The measure- 
ments reviewed here determine diffusion coefficients, not 
the physical nature of intermolecular forces, so our em- 
phasis is on identifying the class of model, not the type 
of force, that is significant for solutions. 

Comparisons were made via non-linear least-squares 
fits. The quantity minimized by the fitting algorithm 
was the mean-square difference between the data and the 
fitting function, expressed as a fraction of the value of 
the fitting function. This quantity is the appropriate 
choice for minimization if the error in the measurement 



7 




c (g/L) 

FIG. 1: D s of 64.2 kDa dextran in water, as obtained with 
PFGNMR by Brown, et al4-, and a stretched-exponential fit. 



is some constant fraction of the value of the quantity 
being measured, e.g., if regardless of the value of D s , D s 
was measured to within 1%. In some cases, one or more 
potentially free parameters were held constant ( "frozen" ) 
during the fitting process. For each fit, the corresponding 
Table reports the final fitting error. 



III. SELF-DIFFUSION 

This Section presents measurements of the true self- 
diffusion coefficient, which describes the motion of a 
labelled chain through a solution of elsewise identical 
chains. Measurements were primarily made with pulsed- 
field gradient nuclear magnetic resonance (PFGNMR) 
and forced Rayleigh scattering (FRS). Results are pre- 
sented alphabetically by first author, showing for each 
paper the data and fits to stretched-exponential forms. 
Fitting parameters appear in Table [I] A more detailed 
analysis of the fitting coefficients appears in Section VI. 
Whenever possible, polymer concentrations have been 
converted to grams/liter. 

Brown, et al4i report self-diffusion and sedimentation 
coefficients for dextran (M n = 44 kDa, M w = 64.2 kDa) 
in water, using PFGNMR to determine D s . They report 
D s and s as functions of c for concentrations as large as 
250 g/L. Figures ^ shows their data and the correspond- 
ing stretched-exponential fit. The self-diffusion data are 
described accurately by the stretched-exponential form. 

Brown and the same collaborators^ used PFGNMR to 
measure D s of of narrow (M w /M n of 1.02-1.20) polyethy- 
lene oxides in water. Polymer molecular weights were in 
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FIG. 2: Self-diffusion coefficients of (top to bottom) 73, 148, 
278, and 661 kDa polyethylene oxides in water, using data of 
Brown, et al^-, and exponential fits. 
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FIG. 3: D a of [top to bottom] 2, 110, 233, and 350 
kDa polystyrenes in CCL as obtained with PFGNMR by 
Callaghan and Pindeu— ■ *— , and fits to stretched exponen- 
tials. 

the range 73-661 kDa; polymer concentrations ranged up 
to 70 g/1. The same paper reports measurements on these 
systems of the mutual diffusion coefficient (from QELSS) 
and the sedimentation coefficient. Figure [21 shows the 
measurements of D s . Because measurements were only 
reported over a limited concentration range, the data 
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FIG. 4: D s of [top to bottom] 110, 233, and 350 
kDa polystyrenes in CeD6 as obtained with PFGNMR by 
Callaghan and Pindeu— *— , and fits to stretched exponen- 
tials. 



FIG. 6: D 3 of 125 kDa polystyrene in toluene as obtained 
with PFGNMR by Fleischej^, and a fit to a stretched expo- 
nential. 
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FIG. 5: D s of [top to bottom] 262, 657, and 861 kDa 
polystyrenes in cyclopentane near the theta temperature, as 
obtained with FRS by Deschamps and Leg and fits to 
stretched exponentials. 

were fit both to a pure [y = 1) and to a stretched ex- 
ponential in c. There is excellent agreement between D s 
and a pure exponential in c, and almost no improvement 
in the quality of the fit attendant to allowing v to be a 
floating parameter in the fit. 

Callaghan and Pinder 43 i 44 i 45 used PFGNMR to study 



D s of 2, 110, 233, and 350 kDa polystyrenes with M w /M n 
in the range 1.06-1.10. Polymer concentrations were as 
large as several hundred g/L. The solvents were CCU 
and hexadeuterobenzene. In CCI4 concentrations well 
into the dilute regime were observed; in GqDq the mea- 
surements of D s were not extended to low concentra- 
tion. The data and corresponding stretched-exponential 
fits appear in figs. and 0] For each polymer:solvent 
pair, a stretched exponential with constant parameters 
describes D s (c) well all the way from the lowest to the 
largest concentrations studied. The prefactor a tends to 
increase with increasing polymer M, while D a and the 
exponent v generally decrease with increasing polymer 
M. 

Deschamps, et al4& used FRS to study self-diffusion of 
polystyrenes in cyclopentane in the vicinity of the theta 
point. The polymer molecular weights were 262, 657, 
and 861 kDa; polymer concentrations ranged from 1 to 
240 g/L. The polymer polydispersity was in the range 
M w /M n « 1.1 — 1.3. Figure [S] shows their data and 
fits to stretched exponentials. For the 657 kDa polymer, 
D s was only reported over a narrow, elevated concen- 
tration range. For the 657 kDa polymer, the stretched- 
exponential fit was therefore made by interpolating D Q 
and v from their values for the 262 and 861 kDa polymers, 
leaving a as the only free parameter. For all three poly- 
mer molecular weights, the stretched-exponential forms 
are in good agreement with experiment. The stretched- 
exponential form for D s (c) remains valid after a change 
from good to near-theta solvent conditions. 

As an aside, the data of Deschamps, et al4P illustrate 
well the principle that an experimental test of a par- 



9 




1 10 100 



c (g/L) 

FIG. 7: D s of [top to bottom] 15, 530, and 730 kDa poly- 
dimethylsiloxane in toluene, as obtained with PFGNMR and 
reported by Giebel, et all— based in part on work of Skirda, 
et al— , and fits to stretched exponentials. 




c (g/L) 

FIG. 8: D s of [from top to bottom] 78, 123, 245, 599, and 745 
kDa polystyrene in benzene, and fits to stretched exponentials 
(fit patameters, Table [I] using data of Hervet, et all— and 
Leger, et al— . 



ticular theoretical model is sometimes not optimal as a 
test for a different theoretical model. Deschamps and 
Leger's objective was to search for deGennes-type scal- 
ing behavior^ 2 of polymer self-diffusion in a theta solvent. 
Scaling behavior is only expected in the semi-dilute con- 
centration regime c > c*. In Deschamps, et al.'s systems, 
the semi-dilute regime was expected to be found only for 
c > 50 — 100 g/L. In the context of their objective of 
studying scaling behavior, there was no strong reason for 
Deschamps, et al. to measure D s at lower concentrations, 
so they rationally did not do so for the 657 kDa system. 
However, an accurate stretched-exponential fit requires 
measurements of D a at small as well as large concen- 
trations. Through no fault of the original authors, the 
range of concentrations studied for the 657 kDa polymer 
restrains the utility of the fits that can be made to some 
of their data. 

Fleischer— used PFGNMR to observe self-diffusion of 
125 kDa polystyrene, M w /M n w 1.02, in toluene for con- 
centrations 80-320 g/L. These concentrations were esti- 
mated to cover 0.5c e < c < 2c e , where c e is the concentra- 
tion above which entanglements were said to be present. 
c e was estimated by applying a deGennes-type model to 
rheological data^ . Fleischer's measurement appear in 
Fig. As seen in the Figure, the measured D s (c) is in 
good agreement with a stretched-exponential form. 

Fleischer's PFGNMR measurements 4 ^ of the incoher- 
ent dynamic structure factor Si nc (q,t) of this system 
show only a single fast relaxation, even under condi- 
tions under which QELSS reveals that a slow mode dom- 
inates the QELSS spectrum. Fleischer's observation that 



the QELSS slow mode cannot be seen in the PFGNMR 
data is immediately reminiscent of results of Zero and 
Ware^ on fluorescence recovery after photobleaching in 
poly-L-lysine solutions at low salt concentration. In each 
set of results, the QELSS spectrum shows distinct slow 
and fast modes. The slow mode becomes dominant un- 
der conditions (higher polymer concentration, lower ionic 
strength) that enhance 'glassy' behavior. In each case, 
the appearance and then dominance of the slow mode, 
when 'glassy' behavior is enhanced, has no effect on the 
self-diffusion of individual chains. This behavior is hard 
to understand if the QELSS slow mode is interpreted as 
arising from long-lived chain clusters with a fixed list of 
members. The behavior is, however, understandable in 
terms of a slow mode arising from the appearance of long- 
lived dynamic structures within which individual chains 
only have short residence times. Such dynamic properties 
of D s and D m for interacting interpenetrating particles 
are also seen in Johnson et al.'saSiai model glasses. 

Giebel, et alM report D s of 15, 530, and 730 kDa poly- 
dimethylsiloxane in toluene as obtained with PFGNMR, 
based in some part on data of Skirda, et al— . The orig- 
inal measurements cover the concentration range 2-900 
g/L of polymer. Fitting parameters are in Table I. Fig- 
ure0shows the actual data and their fits. For each poly- 
mer molecular weight, D s is described well by a stretched 
exponential in concentration. 

Hadgraft, et alM. used QELSS to measure the diffu- 
sion of polystyrenes in benzene as a function of molec- 
ular weight for 24.8 < M < 8870 kDa at 25 C and 
very low polymer concentration. This data is of spe- 
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FIG. 9: D of 245 kDa polystyrene in 245 and 599 kDa 
polystyrene:benzene (open, filled circles, respectively), and 
599 kDa polystyrene in 599 and 1800 kDa polystyrene:benzene 
(open, filled squares), using data of Leger, et alf— and 
patameters in Table Q 



FIG. 10: D s of 130 and 180 g/L polystyrene in dibutylph- 
thalate as a function of molecular weight, showing the 
lower-molecular-weight stretched-exponential and the higher- 
molecular-weight power-law molecular-weight dependences of 
D s , using data of Nemoto, et al.^i, Table II. 



cific interest here in that it supplies values for D s in the 
low-concentration region, as a supplement to data sets in 
which D s was only obtained at relatively elevated con- 
centrations. Omitting the largest-M point (and thereby 
reducing the RMS fractional error in the fit from 18% to 
8.6%), their data follow D Q = 4.54 • lO^M" - 588 . 

Hervet, et al££ and Leger, et al£& studied D s of 
polystyrene in benzene using FRS to measure the relax- 
ation of photoexcitation patterns. Self-diffusion coeffi- 
cients were obtained for polymers of molecular weight 
78.3, 123, 245, 599, and 754 kDa, with M w /M n « 
1.06 — 1.12, for concentrations up to 550 g/L. The same 
technique was used to obtain diffusion coefficients of la- 
belled 245 kDa probe chains in a 599 kDa polystyrene ma- 
trix, and 599 kDa polystyrenes in a 1800 kDa polystyrene 
matrix. Matrix concentrations ranged from 5 to 400 g/L. 
To supplement these measurements, many of which were 
made at elevated polymer concentrations, we used the 
dilute-solution self-diffusion measurements of Hadgraft 54 
to estimate D s at very low c. 

Figure [S] displays the self-diffusion data of Hervet, et 
alm^ and Leger, et ali£& and the fits to stretched expo- 
nentials. Except for the 245 (filled circles) and 598 (open 
squares) kDa polystyrenes, there is good agreement be- 
tween the reported D s (c) and the fits. In the two anoma- 
lous systems, D s at first increases with increasing c and 
then decreases at larger c. The initial increase in D s with 
increasing c appears to be substantially larger than the 
random error in the measurements as inferred from ran- 
dom scatter in the data. The non-monotonic dependence 




c (Volume Fraction) 

FIG. 11: D s of poly(ethylene oxide), [from top to bottom, 
at <f> ~ 0.1] namely 2kDa in CHCI3 and in dioxane, 20 kDa in 
benzene and in dioxane, 40 kDa in chloroform and in dioxane, 
and 3600 kDa in dioxane, after Skirda, et all— Fig. la. Lines 
are fits to stretched exponentials in volume fraction. 
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FIG. 12: D s of polystyrene, namely [from top to bottom] 
240kDa in benzene and 1300 kDa in benzene and carbon tetra- 
chloride, respectively, after Skirda, et all— ■ Fig. lb. Lines are 
fits to stretched exponentials in the volume fraction <f>. 



of D s on c is unique to this specific set of measurements. 
The lack of agreement between the stretched exponential 
form and the measured D s is therefore ascribed to some 
unique feature of this data, and not to a general property 
of polymers in good solvents. 

Figure El presents D p (c) for the two probe-matrix sys- 
tems and compares with D s (c) of polystyrenes of simi- 
lar molecular weights. For both probe-matrix systems, 
D p (c) is described well by a stretched exponential, with 
no indication of the non-monotonic concentration depen- 
dence seen for D s (c) of the 245 and 598 kDa chains. 
The 245 kDa polystyrene diffuses approximately equally 
rapidly through 245 and 599 kDa matrices. The 599kDa 
polystyrene diffuses markedly more slowly through the 
1800 kDa matrix than through the 599 kDa matrix. 

Nemoto, et alSused FRS to measure the self-diffusion 
coefficient of polystyrene as a function of polymer molec- 
ular weight at fixed concentration. The solvent was 
dibutylphthalate. Polymer molecular weights covered a 
range 43.9 < M w < 5480 kDa with polymer concen- 
trations of 130 and 180 g/L. The polymers were quite 
monodisperse, with 1.01 < M w /M n < 1.09, except for 
the 5480 kDa material, for which M w /M n rj 1.15. 

Nemoto, et al.'s results^ appear in Fig. ^] A 
stretched-exponential molecular weight dependence of D s 
is not observed at all M. At the two concentrations 
studied by Nemoto, et al., the molecular weight depen- 
dence of D s has a transition near M ~ 800kDa. At 
lower molecular weights, a stretched exponential in M 
describes D S (M) extremely well. At larger molecular 
weights, one finds a power-law dependence D s ~ A/~ 7 
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FIG. 13: D 3 of hydrogenated polybutadienes in alkanes, with 
polymer molecular weights [from top to bottom] 4.9, 10.3, 
23.3, 53.2, 111, 364, and 440 kDa, and the fitted scaling form 
D s = 55.O0 _1,7 M -2 ' 42 , showing that three parameters and 
the scaling form suffice to describe D s (c, M) over a wide range 
of both parameters. Measurements are from Tao, et all— +. 

with 7 w 2.49 at 130 g/L and 7 » 2.22 at 180 g/L. 

Nemoto, et al^£ used FRS and a cone and plate vis- 
cometer to determine D s and the steady-state shear vis- 
cosity 77 of concentrated solutions (40 and 50 wt%) of 44 
and 355 kDa polystyrene in dibutylphthalate. In terms 
of the transient lattice models, the solutions of the 44kDa 
polymer are expected to be unentangled, while the solu- 
tions of the 355 kDa polymer are expected to be entan- 
gled. They found that D s /T and the fluidity 77 _1 have 
virtually the same dependence on temperature, at either 
concentration, both for the 44 kDa polystyrene and for 
the 355 kDa polystyrene. 

Skirda, et alM used PFGNMR to study the self- 
diffusion of polyethylene oxides (M = 2, 20, 40, and 
3000 kDa) and polystyrenes (M n = 240 and 1300 kDa) 
in chloroform, benzene, dioxane, and carbon tetrachlo- 
ride over a full range of polymer volume fractions <j>. 
M w /M n was w 1.1 for the polyethylene oxides (except 
the 3000 kDa polymer, for which M w /M n ss 2) and w 1.2 
for the polystyrenes. Figures ^2 an d El show results 
for PEO and polystyrene, respectively, and the match- 
ing stretched-exponential fits, using parameters from Ta- 
ble I. For each polymensolvent combination, a stretched- 
exponential form fits the data well. 

Tao, et alM measured D s (using PFGNMR and for- 
ward recoil spectroscopy) and 77 (from the dynamic shear 
moduli) of hydrogenated polybutadienes in alkane sol- 
vents. Polymer volume fractions <f> extended from 0.2 up 
to the melt while polymer molecular weights cover two 
orders of magnitude, from 4.8 to 440kDa, with M w /M n 



12 




0.2 0.4 0.6 0.8 1 
c (Volume Fraction) 

FIG. 14: Same as Fig. 1131 except lines now show the fitted 
form D s = 5.54- 10 3 M -2 ' 42 exp(-5.026c°' 5 M 00019 ). RMS 
fractional errors in the fits shown in Figs. ll3l and ll4l are equal; 
a stretched exponential in c and M thus suffices to describe 
this data. 



of 1.01-1.03 or less. Tao, et al. fit their self-diffusion data 
to a scaling description D s ~ <fi a M b . When they^P forced 
a = —1.8, a one-parameter fit found that the averaged 
Dgcj) 1 - 8 is ~ M~ 2A1 . 

Figures IT31 and ITU show Tao, et al.'s data as fit to to 
power-law and stretched-exponential forms, respectively. 
In Fig. EH all points were fit simultaneously to a scal- 
ing equation, finding D s ~ c M 2A2 , with a fractional 
root-mean square fitting error of 20.2%. These exponents 
differ slightly from those of Tao, et alS: Tao, et al. did 
sequential 2-parameter fits first to determine a and then 
to determine 6, while we did a single three-parameter fit 
to all points. FierureHHshows a fit of all data points to the 
stretched-exponential form D oa M~ z exp(— ad^M 1 ). The 
factor M~ z appears here because we are combining data 
on polymers with multiple molecular weights, and the ex- 
trapolations c — > of measurements at different molec- 
ular weight should extrapolate to a different D a at each 
M. Formally, D OD is the extrapolated zero-concentration 
diffusion coefficient of a highly hypothetical polymer hav- 
ing a molecular weight of unity. On forcing v — 0.5, we 
obtain 7 = 1.9 ■ 10" 4 and z = 2.42 with an RMS frac- 
tional error of 20.2%. Treating yasa free parameter finds 
v w 0.24, with virtually the same value of z, 7 ~ 0, and 
only a slightly improvement (to 19.6%) in the fit error. 
With either value for v, the molecular-weight dependence 
of D s is almost entirely determined by the prefix M~ z . 
The exponential itself has only a negligible dependence 
on M. 

In Tao, et alS's systems, scaling-law and stretched- 
exponential forms for D s thus provide equally good de- 
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FIG. 15: D s of xanthan, molecular weights [from top to 
bottom] 0.45, 0.99, 1.9, 3.8, and 9.4 MDa, in water, and fits of 
the lower-concentration data to stretched exponential forms. 
For the 3.8Mda polymer, solid and dashed lines represent fits 
to the first 7 and 8 data points, respectively. Measurements 
are from Tinland, et all— Fig. 3. 



scriptions of the concentration dependence of D s . The 
scaling and stretched-exponential forms also agree as to 
the molecular weight dependence of D s at fixed c, namely 
D s ~ M~ 2A . Tao, et alS concluded that a scaling-law 
description of their data is correct. The analysis here cor- 
roborates this statement, but shows that it is incomplete, 
in that stretched-exponential forms describe equally ac- 
curately the measured D s (c, M). 

Qualitatively, the larger part of Tao, et al.'s viscos- 
ity measurements^ might have been expected to be in 
the larger-concentration meltlike (77 ~ c x ) rather than 
the smaller-concentration solutionlike (77 ~ exp(ac y )) 
regime. The transition in the concentration dependence 
of rj between these two functional forms is not transpar- 
ently evident in D s (c), whose concentration dependence 
is consistent with a stretched exponential in c for con- 
centrations up to the true melt. However, contrary to 
those other systems reviewed in this article, from which 
molecular-weight dependences can be extracted, here the 
molecular weight exponent 7 of exp(— ac v M 1 ) is very 
nearly zero, so in Tao, et al.'s systems the concentration 
dependence of D s is very nearly the same at all M. 

Tinland, et ali^i report on self-diffusion of xanthans 
of molecular weight 0.45-9.4 MDa, with M w /M n in the 
range 1.2-1.4, at concentrations 0.01-40 g/L. Xanthan 
forms wormlike chains. The experimental data of Tin- 
land, et al. obtained with FRAP thus differs from D s 
of almost all other polymers, in that dln(D s )/d\n(c) 
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FIG. 16: D s of / = 2 (open circles) and / = 8 (closed circles) 
polyisoprenes in CeFsCl, both with M apa n = 5kDa, and fits 
to stretched exponentials in c. Data from von Meerwall et 
alS, Fig. 1. 
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FIG. 17: D s from von Meerwall, et ahS on [from top to 
bottom] / = 2, 3, 8, 18 polyisoprenes with M spa n = 5kDa and 
f — S polyisoprene with M 3pan = 14kDa in CCU, and fits to 
stretched exponentials in c. 



does not decrease monotonically with increasing c. This 
phenomenon, which was noted by the original authors, 
is plausibly related to the appearance- 1 - of a lyotropic 
liquid-crystalline phase in this material at elevated c. 
None of the polymer dynamics models discussed above 
would be expected to remain valid while the system un- 
derwent a phase transition, so the behavior observed by 
Tinland, et al. does not contradict any model for polymer 
dynamics. Figure IT51 shows Tinland, et al.'s** measure- 
ments and fits of the lower-concentration data at each 
molecular weight to a stretched-exponential form. As 
seen from Table I, the somewhat large RMS fractional 
errors show that agreement between the measured data 
and the functional forms is not outstanding. The fit de- 
pends marginally on the number of data points included 
in the analysis. For the 3800 kDa polymer, we indicate 
(solid, dashed lines) the fits to the first 7 or 8 data points. 

For all but the largest-M polymer, there is a con- 
centration c** at which D s (c) deviates from its low- 
concentration decline. In Tinland, et al.'s language^!, c** 
is the concentration of the higher-concentration bound- 
ary of the semidilute regime. From Fig. 1151 c** appears 
to decrease with increasing polymer molecular weight. 

von Meerwall, et al. 62 used PFGNMR to measure D s of 
linear and star polyisoprenes over a near-complete range 
of concentrations (polymer weight fraction 0.01 < x < 1). 
The arm number / ranged from 2 to 18. Molecular 
weights of single arms were 5 and 18 kDa. Figure ^] 
compares von Meerwall, et al.'s§£ data on / = 2 and 
/ = 8 polyisoprene stars having 5kDa arms. Figure 
1171 shows D s for / = 2,3,8, and 18 star polyisoprenes 



in CCI4. In the original paper, data were reported as 
smooth curves, not as measured points. To create Fig. 
1171 the smooth curves were sampled; fits were made to 
the sampled points. From Fig. El a t &H concentrations 
increasing the arm count at fixed span molecular weight 
reduces D s : Increasing four- fold the number of arms re- 
duces D s by a factor of 2 to 3. At / = 8, an increase in 
span molecular weight reduces D s . 

In both figures, the solid lines are fits to stretched ex- 
ponentials. Without exception, D s (c) for each / and M is 
described well by a stretched exponential in c. From the 
fit parameters in Table I, for fixed M span and increasing 
/ one finds that D Q and v fall while a increases, a two- 
fold decrease in D Q via increasing / being accompanied 
by a 20-fold increase in a. For the smallest molecular 
weights (10-16 kDa) studied, v > 1 is observed. 

von Meerwall, et aliSi used PFGNMR to study lin- 
ear and 3-armed star polybutadienes and polystyrenes 
in CCI4. They report low-concentration data on all sys- 
tems, and extensive concentration dependence measure- 
ments on some star polymers. Polymer molecular weights 
ranged from 2.3 to 281 kDa with polydispersities M w /M n 
of 1.03-1.07. As seen in Figs. 1181 and !19l for most systems 
the observed concentration range afforded a one-order-of- 
magnitude variation in D s . 

Figure |20 gives Ref. l63T s measurements on the 6.5, 
8.3, 29, and 76 kDa 3-armed stars, which were made 
over a far wider concentration range that afforded a 2.5- 
order-of-magnitude variation in D s . Fits of these data to 
stretched-exponential forms describe well D s (c). As ini- 
tially noted by von Meerwall, et ali££, 'the slopes [in Fig. 
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FIG. 18: D s of [from top to bottom] 2.3 kDa linear, and 6.5, 
16.1, 21, and 26 kDa three-armed polybutadienes in CCk, and 
fits to stretched exponentials, using data of von Meerwall, et 
alS, Fig. 2. 



|2*U] change continuously' '. . . deGennes' prediction of a 
concentration scaling regime D s ~ c~ (c* < c < c**) 
is not borne out by our data ... at any molecular weight'. 

von Meerwall, et al. 64 used PFGNMR to measure self- 
diffusion of 10, 37.4, 179, 498, and 1050 kDa polystyrene 
in tetrahydrofuran at concentrations 6-700 g/L. The 
same technique was also used to measure D s of tetrahy- 
drofuran and hexafiuorobenzene in the same polymer so- 
lutions. Von Meerwall, et al. did not report D 8 for their 
polymers in dilute solution. D s for these systems in the 
dilute limit is therefore inferred here from the molecular 
weight dependence of D s observed by Hadgraft, et al^ 
for polystyrene in benzene, together with the viscosities 
of tetrahydrofuran and benzene, von Meerwall, et al.'s 
measurements, and stretched-exponential fits appear in 
Fig. |^ using parameters given in Table I. 

Wesson, et al.® 5 used FRS to measure self-diffusion 
of polystyrenes in tetrahydrofuran and benzene. The 
polystyrenes had M of 32, 46, 105, 130, and 360 kDa, and 
were observed for concentrations in the range 40-500 g/L. 
Wesson, et al.'s measurements were here supplemented 
by extreme low-concentration points calculated from re- 
sults of Hadgraft, et ali 5 ^. These points were included 
in the fits as having been taken at g/L. Experimental 
results and the corresponding stretched-exponential fits 
appear in Fig. |2H Because D s here covers four orders of 
magnitude, on the scale of the figure the fits look very 
good. In fact, RMS fractional errors are in the range 15- 
21%, making these among the poorer fits in this section 
to a stretched exponential. 

Xuexin, et al^£ used PFGNR to measure D s of lin- 




25 50 

c (g/L) 

FIG. 19: D 3 of [from top to bottom] 75, 76, 90, 161, 227, and 
281 kDa polybutadienes in CCk, the 75 and 90 kDa polymers 
being linear and the others being three-armed stars, and fits 
to pure exponentials, using data of von Meerwall, et alp—, Fig. 
3. 




c (g/L) 

FIG. 20: D a of [from top to bottom] 6.5 8.7, 29, and 76 
kDa three-armed polybutadienes in CCk, and fits to stretched 
exponentials, using data of von Meerwall, et al^, Fig. 4. 

ear and 18-armed star polyisoprenes in CCk over a wide 
range of c and a 100- fold range of M. Their results (with 
concentrations replotted in g/L) appear as Figs. [23 EH 
and 1251 In every case the concentration dependence of 
D s is described well by the stretched exponential form, 
with fitting parameters given in Table I. 
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FIG. 21: D s of 10, 37, 198, 498, and 1050 kDa linear 
polystyrenes in tetrahydrofuran, and fits to stretched expo- 
nentials in c. Data after von Meerwall et alp—, with supple- 
mental low-concentration data based on Hadgraft, et alp—. 



FIG. 22: D a of [from top to bottom] 32, 46, 105, 130, and 360 
kDa polystyrene in tetrahydrofuran, and associated stretched 
exponential fits. Data are from Wesson, et alp—, Table II and 
associated Figures, as supplemented by the low-concentration 
measurements of Hadgraft, et alf*— . 



Figure [23] shows D s (c) for four low-molecular- weight 
(61, 92, 193, 216 kDa) 18-arm polyisoprenes. With in- 
creasing M and fixed /, D a and v decrease, while a 
increases, the increase in a being modestly more rapid 
with increasing M than is the decrease in D . These 
trends are continued in Fig. [2J] which shows data from 
ref.Eifor larger-M (344, 800, 6300 kDa) stars. The dis- 
played data on 302 kDa linear polyisoprene are fit by very 
nearly the same D a , a and v as is the data on a much 
larger (800 kDa) 18-arm star polyisoprene. Finally, Fig. 
1251 shows Xuexin, et al.'s data on D s (c) for three linear 
polyisoprenes, molecular weights of 70.8, 251, and 302 
kDa, albeit over a narrower range of c than in Figure 1231 

Nine of the above papers, namely Brown and Stilbs^, 
Callaghan and Pinder— , Deschamps and Leger—, Giebel, 
et al—, Leger, et al—, von Meerwall, et al— , Wesson, 
et al— , and Xuexin, et al— report self-diffusion coeffi- 
cients at a series of concentrations and polymer molecular 
weights for a series of homologous polymers. A simulta- 
neous fit of each of these sets of measurements to a joint 
function of c and M is then practicable. Here fits were 
made to eq. 1171 using the convention that the molecular 
weight of the sole polymer species is P. The fits forced 
5 = to eliminate a notional dependence of D s on the 
molecular weight M of a non-existent second polymer. 
Results of the fits appear in Table IV and Figs. 1261251 

In a majority of cases, the fits were quite good, with 
RMS fractional errors in the range 6-18%. Fits to the 
measurements of Leger, et al— and Wesson, et al— were 
less satisfactory; these are discussed separately below. 
Xuexin, et al— cover an extremely broad range of M in 



their study of 18-armed stars. Over this range, v from 
fits to data sets covering chains with a single molecular 
weight changes substantially, so a fit of all data over a full 
range of M to eq. El with v a fixed constant, does not 
work well for very small or very large molecular weights. 
On limiting the fit to intermediate values for P, good 
agreement between the data and the fitted forms are en- 
countered. Solid curves in Fig. 1271 were calculated using 
the parameters to the fit to data in the intermediate-P 
range. 

Omitting momentarily the fits to results of Leger, et 
al— and Wesson, et al— : The exponent a for the molec- 
ular weight dependence of the bare diffusion coefficient 
is consistently -0.5. Except for Browne, et al.'s work on 
polyethylene oxide: water, the concentration exponent 
v is in the range 0.5-0.75; Brown, et al.'s data^ imply 
v w 0.93. The molecular weight exponent 7 is in the 
range 0.32-0.46, again with the exception of fits to data 
of Brown, et al— , for which 7 sa 0.6. Inspection of the 
figures indicates that a joint stretched exponential in c 
and M fits each data set well, with no systematic devia- 
tions for particular values of c or M. 

The two data sets that are fit less well by eq.H7lappear 
in Figure 1251 The merged fit to Leger, et ally's data on 
polystyrene:CCi4 is poor. However, Leger, et al's data 
shows features - notably a non-monotonic dependence 
of D s on c - that appears in data on no other polymer 
system, including other experiments that determined the 
concentration dependence of D s of the same polymer. 
We infer that the poor fit of eq.Elto Leger, et al.'s data 
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FIG. 23: D a of [from top to bottom] 61, 92, 193, and 216 
kDa 18-arm star polyisoprenes in CCk, and fits to stretched 
exponentials, using data of Xuexin, et all— , Fig. 3. 



FIG. 25: D s of [from top to bottom] 70.8, 251, and 302 kDa 
linear polyisoprenes in CCk, and fits to stretched exponen- 
tials, using data of Xuexin, et al^, Fig. 4. 




FIG. 24: D s of 302 kDa linear polyisoprene (filled circles) and 
[from top to bottom] 340, 800, and 6300 kDa 18-armed star 
polyisoprenes (all in CCk) and fits to stretched exponentials, 
using data of Xuexin, et all— , Fig. 2. 



arises from features unique to this polymer sample and 
set of measurements and not to a generic behavior of D s 
for polystyrene solutions. 

The merged fit to Wesson, et al.'s^ measurements is 
also poor. However, this set of data is limited to elevated 
concentrations in which D s /D s q < 1, generally substan- 



tially without matching low-concentration measurements 
on the regime in which D s « D s q. By inspection of Fig. 
128b . the merged fit works reasonably well at small and 
large polymer molecular weight, but is least satisfactory 
at intermediate molecular weights. In particular, the 105 
and 130 kDa polymers have nearly the same P but at ele- 
vated concentrations substantially different self-diffusion 
coefficients. 



IV. DIFFUSION OF PROBE CHAINS 
THROUGH MATRIX POLYMER SOLUTIONS 

This Section reviews measurements on the diffusion 
of polymeric probe molecules through solutions of a dif- 
ferent polymer. These experiments involve intrinsically 
ternary solutions in which the molecular weight P of 
the probe polymer and the molecular weight M of the 
matrix polymer are not the same. In some cases, the 
probe and matrix polymers in the solution have a com- 
mon monomer, and differ only in their molecular weights. 
In other cases, the probe and matrix polymers are chem- 
ically distinct. Studies are again presented alphabeti- 
cally by first author, together with fits of the data sets to 
stretched exponentials in concentration and chain molec- 
ular weights. Fitting parameters appear in Tables HTI and 

nn 

Brown and Rymden§£ used quasielastic light scatter- 
ing to study the diffusion of linear polystyrenes and 
coated silica spheres through polymethylmethacrylate in 
toluene. Toluene and PMMA are almost exactly index- 
matched, so scattering from these systems was dominated 
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FIG. 26: D s as measured by (a) Brown, et al.— (cf. Fig.0, (b) Callaghan, et al.— (cf. Fig. El, (c) Callaghan, et al4* (cf. Fi; 
^J, and (d) Deschamps, et al4- (cf. Fig.QjJ, and fits to eg. 1171 leading to the parameters in Table ITTTl 



by scattering from the dilute probe chains. The matrix 
PMMA's had molecular weights in the range 110 kDa- 
1.43 MDa. Probe polystyrenes had molecular weights of 
2.95, 8, and 15 MDa, with M w /M n of 1.06, 1.08, and 1.30, 
respectively. In the original paper, data were reported af- 
ter normalization by an unspecified diffusion coefficient 
obtained in the absence of the matrix polymer. 

Reference |H3 reported how D p depends on matrix 
concentration and molecular weight. Figure I2$k shows 
Dp/Dpo for the 8MDa polystyrene diffusing through each 
of six matrix polymethylmethacrylates. Each solid line 
represents a fit to a stretched exponential in matrix poly- 
mer concentration, using parameters in Table II. For 



Dp/Dpo < 10~ 3 or so, a condition attained only with 
the two largest matrix polymers, the measured D p devi- 
ates markedly downward from a stretched exponential, so 
that fits to all data on these two systems show fractional 
RMS errors much worse (32-45% rather than 5-8%) than 
fits to the same probe polymer with the smaller matrix 
polymers. On excluding the final few points from the two 
fits, curves with far lower RMS fractional errors were gen- 
erated. These are the curves shown in the Figure. It is 
difficult to determine from these measurements whether 
there is a systematic change from stretched-exponential 
to some other concentration dependence at very small 
Dp/DpQ, or whether the apparent deviations arise from 
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FIG. 27: D s as measured by (a) Giebel, et al&, (b) von Meerwall, et al?S and (c), (d) Xuexin, et al.^, and fits to eg. IT71 
leading to the parameters in Table HTT1 Other details as in Figs. andlH respectively. 



experimental challenges at very small D p . 

Figure 129b shows the same data, now fit simultane- 
ously to eqn. El the joint stretched exponential in c and 
M. Fitting parameters appear in Table IV. Except for 
the lowest- M matrix polymer, the fits are almost as good 
as the individual fits shown in Figurel2*9"k. For the lowest- 
M 101 kDa matrix polymer, the fitted form predicts too 
strong a dependence of D p upon c. 

Brown and RymdenSZ, also examined how D p depends 
on probe molecular weight. Figure QUI shows D p of the 3, 
8, and 15 MDa probe polystyrenes, all diffusing through 
the 445kDa PMMA. Fits are to stretched exponentials 
in c, leading to parameters given in Table [H] The ma- 



trix chains are all much larger than the probe chains. As 
noted by the original authors^, the three curves come 
very close to superposing except perhaps at the very high- 
est matrix concentrations examined. 

Brown and Stilbs^ used PFGNMR to measure the 
probe diffusion coefficient of polyethylene oxide in aque- 
ous solutions of dextran. Polyethylene oxides had molec- 
ular weights of 73, 278, and 1200 kDa with M w /M n of 
1.02-1.12; dextrans had molecular weights of 19, 110, and 
510 kDa. The 1200 kDa PEO represented the lower limit 
at which D p could be determined with then-available 
technology; the authors limited their detailed analysis 
to the two lower-molecular-weight probes. 



19 




io u m 



10"' 



10 2 \m 



10"° 



10" 



10 



10 




10 100 

c (g/L) 



-6 

0.001 0.01 0.1 1 

c (g/L) 



10 100 



10" D s 



10"' = 



^2 10" 8 
E 

Q M 10" 9 b 



10" 10 



10 



11 








10 



100 

c (g/L) 



1000 



10 



10"' 



10" 



10 3 



1<T = 



10~ 5 



10" 




0.001 0.01 0.1 1 

c (g/L) 



10 100 



FIG. 28: Joint fits to all data of (a) Leger, et alM and (b) 
Wesson, et all— , showing the poor quality of fits of eg. 1171 to 
these data sets. Fit parameters are in Table ITTTl all other plot 
properties are the same as Figs. 151 and 1221 respectively. 



Figures 131b and 131b show D p of the 73 and 278 

kDa polyethylene oxides, normalized by the measured 
D q of the same probes in the absence of the matrix 
polymers. All measurements were simultaneously fit to 
Dp /D o exp(— ac" P 1 M s ), yielding the smooth curves 
shown in the Figures and the parameters listed in Ta- 
ble IIVI Agreement between the data and the fitted 
curves was very good for the 73 kDa probe. For the 
278 kDa probe in the 19 kDa matrix, the fitting func- 
tion significantly underpredicts D p . From the fitting pa- 
rameters, Dp/ D q has a very weak dependence on the 



FIG. 29: Dp/Dpo of 8 MDa polystyrene diffusing through 
[top to bottom] 101, 163, 268, 445, 697, and 1426 kDa poly- 
methylmethacrylate in toluene—, and (a) separate fits of the 
data for each matrix polymer to a stretched exponential in 
c, and (b) simultaneous fit of all data to a single stretched 
exponential in c and M. 



probe molecular weight (other than that hidden in D q), 
but has a marked dependence on the matrix molecular 
weight. The same data were also fit, individually for 
each probe:matrix pair, to a stretched and a pure expo- 
nential in c. Without exception, for each probe:matrix 
pair D s (c) follows accurately the exponential form, RMS 
fractional errors being in the range 1-4%. Because D p 
varied over a limited range, the parameters reported in 
Table II reflect fits to the pure exponential. 
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FIG. 30: Dp/D p o against c for 3 (open circles), 8 (closed 
circles), and 15 (squares) MDa polystyrenes diffusing through 
445 kDa polymethylmethacrylate in toluene, using data of 
Brown and Rymden— , and fits to stretched exponentials in 
matrix concentration. 



Daivis, et used quasielastic light scattering spec- 
troscopy to measure diffusion of relatively dilute 864 kDa 
dextran in not-necessarily-dilute solutions of 20.4 kDa 
dextran. Polymer polydispersities were in the range 1.24- 
1.3. The concentration of the lower-molecular- weight 
dextran ranged up to 166 g/L. Analysis of the bimodal 
QELSS spectra of these systems shows that the slower 
mode corresponds to probe diffusion by the 864 kDa dex- 
trans. As seen in Figure 02 the data fit reasonably well 
to a simple exponential concentration dependence, using 
parameters found in Table II. 

In a separate paper, Daivis, et ali 7 ^ used QELSS 
and PFGNMR to measure the diffusion of a 110 kDa 
polystyrene, M w /M n = 1.06, through solutions of 
110 kDa polyvinylmethylether, M w /M n 1.3, in the 
PVME's isorefractive solvent toluene. Good agreement 
was found between the QELSS measurements of D p and 
the earlier measurements of D p by MartinSi on the same 
system. 

Figure 13^1 shows Daivis, et al.'s dataZS as obtained us- 
ing both physical methods. In the figure, lines represent 
separate fits of each data set to a stretched exponential 
in matrix polymer concentration, with fitting parameters 
given in Table [H] Stretched exponentials in c describe 
well each data set. The QELSS data is significantly less 
scattered than the PFGNMR data, so the parameters 
from the former's fit are probably to be preferred. 

De Smedt, et al^ use d FRAP to measure the diffusion 
of 71, 148, and 487 kDa dextrans (M w /M n < 1.35), la- 
beled with fluorescein isothiocyanate, through solutions 
of hyaluronic acid (M„ = 390 kDa; M w = 680 kDa). 
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FIG. 31: D p /D o of (a) 73 kDa polyethylene oxide and (b) 
278 kDa polyethylene oxide in [top to bottom] 19, 110, and 
510 kDa dextrans in aqueous solution, and simultaneous fit 
of all measurements to a stretched exponential in c, M, and 
P, using data of Brown and Stilbs— . 



Hyaluronic acid concentrations ranged from the dilute 
up to 18 g/L. D p of the dextrans varied roughly five-fold 
over this concentration range. De Smedt, et al.'s data 72 
appear in Fig. 02] together with stretched-exponential 
fits using the parameters in Table II. As seen in the Fig- 
ure and initially reported by the original authors, the 
data fits well to stretched exponential forms. 

Hadgraft, et alM used QELSS to study the diffu- 
sion of polystyrene probe polymers, molecular weights 
25, 162, 410, 1110, and 4600 kDa through solutions of 
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FIG. 32: D p of 864 kDa dextran in solutions of 20 kDa dex- 
trans and fits to stretched exponentials, using data of Daivis, 
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FIG. 33: D p of 110 kDa polystyrene through 110 kDa 
polyvinylmethylethentoluene, based on QELSS (circles) and 
PFGNMR (squares) measurements of Daivis, et all— ■ , to- 
gether with fits (solid, dashed lines, respectively) to stretched 
exponentials in the matrix concentration. 
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FIG. 34: Dp of [from top to bottom] 71, 148, and 487 kDa 
dextrans in M w 680 kDa hyaluronic acid, and fits to stretched 
exponentials, using data of De Smedt, et air—. 
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FIG. 35: Dp of [from top to bottom] 25, 162, 410, 1110, and 
4600 kDa polystyrenes in 1.05 MDa polymethylmethacrylate: 
CCU as a function of polymethylmethacrylate concentration, 
and fits to stretched exponentials, using data of Hadgraft, et 



105 kDa polymethylmethacrylate in its isorefractive sol- 
vent benzene at PMMA concentrations up to 100 g/L. 
Polystyrene and PMMA are not compatible, implying 
that the radius of the polystyrene chains may have de- 
pended very strongly on the matrix polymer concentra- 
tion. As seen in Fig. [213 D p of the 25 kDa polystyrene was 



substantially independent of PMMA concentration. The 
data on the 410 kDa polymer is significantly more scat- 
tered than is data on the other polystyrenes. The range 
of variation of D p is sufficiently small (roughly a factor 
of three) that the fits to these data are less reliable than 
are the fits to data on some other systems. Nonetheless, 
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ley, et al.'s dataS are given in Fig. l3"fJk . together with 
fits of separate stretched-exponentials in c to each data 
set, yielding parameters shown in Table [H] The entire 
data set was also fit simultaneously to the joint stretched 
exponential of eq. If 71 yielding parameters seen in Table 
11111 and fitted curves seen in Fig. 136b . The lack of mea- 
surements at very low matrix concentration substantially 
broadens the range of fitting parameters that yield rea- 
sonable descriptions of this data. The fit to the joint 
stretched exponential is substantially less satisfactory 
than the individual fits at each P to separate stretched 
exponentials (RMS fractional errors of 34% rather than 
9-24%). The inadequacy of eq. El relates primarily to 
the f 050 kDa probe chains, for which D p was determined 
only over a limited range. As seen in TableHni excluding 
data on the f050 kDa probes from the simultaneous fit 
leads to a marked reduction in the RMS fractional fit er- 
ror (to 24%), modest changes in a and v, but only small 
changes in the other fitting parameters. 

Kent, et ali^ applied static and quasielastic light scat- 
tering to measure radii of gyration and diffusion coeffi- 
cients of 233 and 930 kDa polystyrenes in 7, 66, 70, 840, 
and 1300 kDa polymethylmethacrylates in ethyl benzoate 
(for static light scattering) and toluene (for quasielastic 
light scattering). With one exception (66 kDa PMMA), 
M w /M n was always < 1. f 0. Different probe:matrix com- 
binations were used for static and quasielastic light scat- 
tering. D of the polystyrene was measured as a function 
of polystyrene concentration, and linearly extrapolated 
to the dilute- in-polystyrene limit, thereby obtaining both 
the probe diffusion coefficient D p of the polystyrene and 
the initial linear dependence of D p on polystyrene con- 
centration. 

Figure 03 shows Kent, et al.'si 4 measurements of the 
radius of gyration of 930 kDa polystyrene in solutions 
of 7, 70, or f300 kDa PMMA as a function of PMMA 
concentration, together with fits of R g to a stretched ex- 
ponential 



FIG. 36: D p of [from top to bottom] 50, 179, 1050, and 1800 
kDa polystyrenes in orthofluorotoluene solutions of 60kDa 
polyvinylmethylether against matrix polymer concentration, 
using data of Hanley, et all— and (a) separate fits at each 
P and (b) simultaneous fits at all P, with fitting parameters 
from Tables HT1 and HTT1 respectively. 



Dp of the higher-molecular-weight polystyrenes shows a 
stretched-exponential dependence on PMMA concentra- 
tion, with parameters seen in Table II. 

Hanley, et alS used light scattering spectroscopy to 
examine the diffusion of polystyrenes through the ma- 
trix polymer polyvinylmethylether in its isoreff active sol- 
vent orthofluorotoluene. The polystyrenes had molecular 
weights of 50, 179, 1050, and 1800 kDa. The polyvinyl- 
methylether had M w « 60 kDa with M w /M n w 3. Han- 



R g = RgQ exp(— ac v ) 



(18) 



in matrix concentration c. The fits are good through- 
out, using parameters in Table Ivl Unlike the scaling-law 



prediction R„ 



-0.25 



the stretched exponential form 



shows acceptable behavior down to zero matrix concen- 
tration. 

Kent, et al also measured D p of 233 kDa polystyrene 
in solutions of 66 and 840 kDa PMMA, and 930 kDa 
polystyrene through 840 kDa PMMA. Their experimen- 
tal data is shown in Fig. 138b . together with fits to expo- 
nentials using parameters given in Table HT1 Within ex- 
perimental error, the simple exponential fits with v = 1 
forced are as good as the stretched-exponential fits to the 
data: the former are in the figure. Experimentally, the 
scaling prefactor a depends strongly on matrix molecular 
weight (a 12-fold change in M leads to a two-fold change 
in a) but at most weakly on probe molecular weight. Fig- 
ure EHt shows a fit of the same data to a joint stretched 
exponential in c, P, and M, based on parameters in Table 
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FIG. 37: R g of 930 kDa polystyrene in [from top to bottom] 
7, 70, and 1300 kDa polymethylmethacrylate: ethylbenzoate 
as a function of polymethylmethacrylate concentration, and 
fits to stretched exponentials, using data of Kent, et aim—. 



IIIII The joint fit clearly works well for all three polymer 
pairs. 

Kim, et aliS measured the diffusion of dye-labeled 
polystyrenes through matrix solutions of unlabeled 
polystyrenes in toluene. The objective was to test the 
prediction of some scaling models that D p becomes inde- 
pendent of matrix molecular weight if M/P > 1. The 
probe polystyrenes had molecular weights from 10 to 
1800 kDa; matrix chains had molecular weights from 51 
to 8400 kDa. Ref . [fH also reported the dependence of D p 
on the matrix molecular weight for several probes (51, 
390, 900 kDa) at multiple matrix concentrations for ma- 
trix molecular weights in the range 35-8400 kDa. Poly- 
mer polydispersities were largely < 1.06, with a maxi- 
mum of 1.17. Kim, et alii* also report limited data using 
methyl red as a low-molecular weight probe. As seen in 
Fig. 139b . the authors found that D p becomes substan- 
tially independent of M only if M/P > 3. 

With respect to the models discussed in Section II, the 
published derivations of the stretched-exponential form 
refer to polymer chains whose motions are adequately 
approximated by whole-chain translation and rotation. 
These approximations are only likely to be adequate if 
the probe and matrix chains are of similar size, because 
otherwise the whole-chain motions of probe or matrix 
would effectively sample some of the internal modes of 
the other chain species, whether matrix or probe. At 
fixed c, from the hydrodynamic model D p would have 
a stretched-exponential dependence on M if M/P « 1, 
but might well not have a stretched-exponential matrix 
molecular weight dependence if M/P > 1 or Af/Pc 1. 
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FIG. 38: Dp of polystyrene through polymethylmethacrylate: 
toluene with molecular weight combinations P : M [from top 
to bottom] 233:66, 233:840, and 930:840 kDa as functions of 
polymethylmethacrylate concentration, using data of Kent, 
et ali—, and (a) a separate fit to a stretched exponential for 
each combination and (b)a single fit to a joint stretched ex- 
ponential in c, P, and M. 



Kim, et al.'aS results thus are not necessarily inconsis- 
tent with exponential- type models of polymer dynamics. 

Kim, et al.'s dataS for the domain M/P < 3, as seen 
in Figs.EH were fit to the form D P- a cxp (-ac" P 7 M s ), 
yielding parameters in Table IrTTI As seen in Fig-OHl Dp is 
indeed described well by the joint stretched exponential 
over a wide range of c, P, and M, except in the regime 
M/P > 3 in which the form is not necessarily expected 
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to apply. 

Kim, et al. also studied the concentration depen- 
dence of D s and D p at large M/P. Figure 0U| shows 
Kim, et al.'aS determination of self-diffusion of 900 kDa 
polystyrene in toluene, and our stretched-exponential fit 
to the data, yielding the parameters in Table I. Figurel4*Tl 
shows measurements of probe diffusion, using as probes 
methyl red and 10, 35, 100, 390, 900, and 1800 kDa 
polystyrenes, all through matrix polystyrenes in toluene. 
The molecular weight of the matrix polystyrene ranges 
from 51 to 8400 kDa, varying from data point to data 
point, with M/P > 3 and generally M/P > 6. This 
ratio of M/P was chosen by Kim, et al. based on their 
interpretation of their data that showed D p to be inde- 
pendent of M for M/P > 3. As seen in Figs. EDI and EH 
stretched-exponential functional forms using parameters 
in Table II describe the matrix-polymer concentration 
dependences of D s (c) and D p (c) well over a wide range 
of polymer concentrations and probe molecular weights, 
even though the ranges of P and M are very wide. The 
stretched exponential in c continues to describe well the 
concentration dependence of D p in the large M/P range 
which cq. El does not represent well the P and M de- 
pendences of D s . 

Lodge and collaborators have reported an extensive se- 
ries of studies of probe diffusion in polymer solutions, 
using quasielastic light scattering to measure D p of a 
dilute probe polymer, generally polystyrene, through 
the isorefractive matrix polymer: solvent pair polyvinyl- 
methylether: orthofluorotoluene. Variables studied in- 
clude the probe and matrix molecular weights, the ma- 
trix concentration, and the topology (linear and star) of 
the probe polymers. 

An early letter— of Lodge reports D p of 179 kDa 
and 1.05 MDa polystyrenes through a 50kDa polyvinyl- 
methylether. The data, and corresponding stretched- 
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FIG. 39: D a of polystyrene in matrix polystyrene:toluene 
solutions as functions of matrix molecular weight at various 
matrix concentrations, based on data of Kim, et alf— , Fig. 
3. Probe molecular weights were (a) 51 kDa, (b) 390 kDa, 
and (c) 900 kDa. Matrix concentrations [top to bottom] 
were (a) 10, 20, 50, and 100 g/L; (b) 20, 50, and 100 g/L; 
and (c) 10, 20, 40, and 80 g/L. A single stretched exponen- 
tial D P~ a exp(— ac v P s M 1 ) with constant parameters was 
fit (solid lines) to all data in the figures having M/P < 3. 



exponential fits, appear in Fig. 021 Fit parameters are 
in Table ^ A stretched exponential in concentration 
describes well both data sets. 

Lodge and Wheeled compared the diffusion of lin- 
ear and 3-armed star polystyrenes through a high molec- 
ular weight polyvinylmethylether. Polystyrene molecu- 
lar weights were 422 and 1050 kDa for the linear chains 
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FIG. 40: D s of 900 kDa in toluene, based on data of Kim, et 
al.— -, Fig. 1, and a fit to a stretched exponential. 
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FIG. 41: Probe diffusion of (top to bottom) methyl red, and 
10, 35, 100, 390, 900, and 1800 kDa D s polystyrenes through 
high-molecular-weight (MjP > 3) polystyrenes in toluene, 
based on data of Kim, et all—, Table HTl and fits to stretched 
exponentials. The left axis shows zero-matrix-concentration 
data. 

and 379 and 1190 kDa for the stars; the polystyrenes 
were identified as being 'relatively monodisperse'. The 
PVME had an M w of 1.3 MDa, with an estimated 
M w /Mm ps 1.3. D p was obtained over 1 < c < 100 
g/L in PVME concentration. 

Figure 143b shows Lodge and Wheeler's 
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FIG. 42: D p of 179 kDa and 1.05 MDa polystyrene through 
50 kDa polyvinylmethylether in orthofluorotoluene, based on 
data of Lodged, and fits to stretched exponentials in c. 



measurements^, together with fits of data on each 
probe polymer to a stretched exponential in c. D p varies 
over nearly four orders of magnitude. Over the full 
range, agreement of the data with the functional form is 
very good, with RMS fractional errors of 2.6-10% and 
parameters as seen in Table [H] Figure 143b shows the 
same data, with the linear chains and star polymers 
separately fit to a stretched exponential in c and P. 
As seen in Table IIIII the fractional errors in these fits 
are very nearly as good as the fits to individual probe 
species. The P~ a scaling of the zero-concentration 
diffusion coefficient and the P 1 scaling of the scaling 
prefactor a account for the dependence of D p of linear 
and star polymers on probe molecular weight. 

Lodge and MarklandS used light scattering spec- 
troscopy to measure the single-particle diffusion coeffi- 
cients of tracer 12-armed star polystyrenes through so- 
lutions of 140 kDa polyvinylmethylether, M w /M n w 
1.6, in the isorefractive solvent orthofluorotoluene. The 
polystyrenes had M w of 55, 467, 1110, and 1690 kDa, 
with M w /M n < 1.10. Lodge and Markland estimate for 
the matrix that c* w 20 g/L and c e w 100 g/L. Fig- 
ure El shows Lodge and Markland's data. In the two 
graphs, the solid lines represent, respectively, fits to in- 
dividual stretched exponentials in c and to fits to a joint 
stretched exponential in c and P. Fitting parameters ap- 
pear in Tables |n] and IIIII Parameters in Table [H] differ 
modestly from Table HJ of Ref. El The sets of fits are 
excellent, with RMS fractional individual fits here are ex- 
cellent, with errors of 2-4%; the RMS fractional error for 
the joint fit was 12%. The most notable deviation for 
the joint exponential is for the smallest probe at large 



26 




1 10 100 



c (g/L) 

FIG. 43: D p of polystyrenes [from top to bottom: 379 kDa 
/ = 3 star, 422 kDa linear chain, 1050 kDa linear chain, and 
1190 kDa / = 3 star] through 1300 kDa polyvinylmethylether 
in orthofluorotoluene, based on data of Lodge 77 and fits of a 
stretched exponential to (a) data on each probe separately 
and (b) all data on linear chains (solid lines) and separately 
all data on 3-armed star polymers (dashed lines). 
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FIG. 44: D p of 12-arm star polystyrenes (from top to bottom, 
M w of 55, 467, 1110, and 1690 kDa) diffusing through so- 
lutions of 140 kDa polyvinylmethylether:ortho-fluorotoluene 
based on data of Lodge, et all— , their Table 1. Solid lines 
are fits of (a) a separate stretched-exponential for each probe, 
with parameters in [H] and (b) fits of a single stretched ex- 
ponential in c and P to all data with parameters in Table 

urn 



concentration, where the fitting function underestimates 
D p . 

Lodge, Markland, and Wheeler— used light scatter- 
ing spectroscopy to measure the diffusion of 3-armcd 
and 12-armed star polystyrenes through solutions of 
polyvinylmethylether in its isorefractive solvent orthoflu- 
orotoluene. The 3-armed stars had M w of 379 and 1190 
kDa; the 12-armed stars had M w of 55, 467, 1110, and 



1690 kDa. Polyvinylmethyl ethers used as matrices had 
M w of 140, 630, and 1300 kDa. Polystyrenes all had 
M w /M n < 1.1; the matrix polymers had M w /M n rs 1.6. 
Light scattering measurements were also made of the 
radii of gyration of linear, 3-arm and 12-arm stars with 
molecular weights above lMDa in the presence of 250 
kDa polyvinylmethylether at concentrations as large as 
50 g/L. Lodge, et alJS, reported that they modified their 
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FIG. 45: D p of (a) 379 kDa and (b) 1.19 MDa 3-arm star 
polystyrenes in (from top to bottom, 140, 379, and 1300 kDa) 
polyvinylmethylethers in orthofiuorotoluene, based on data of 
Lodge, et ali—». Solid lines are stretched-exponential fits with 
parameters in Table ITT1 



reported D p to take account of the concentration depen- 
dence of a local friction factor, using a process described 
by Wheeler and Lodged (see below). This modification 
factor was removed from Lodge, et al's^ data, before 
making the analysis here, in order that this data be made 
more strictly comparable with the remainder of the lit- 
erature. 

Figures R31 and l4"6l show Lodge, et al. 'admeasurements 
of Dp for 3-armed and 12-armed stars. As is apparent 
from the figures, for every M : P combination a stretched 
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FIG. 46: Dp of 12-arm star polystyrenes (from top to bottom, 
M w of 55, 467, 1110, and 1690 kDa) diffusing through so- 
lutions of 1300 kDa polyvinylmethylethenortho-fluorotoluene 
based on data of Lodge, et alS. Solid lines are stretched- 
exponential fits with parameters in Table II. 



exponential gives an excellent description of the the con- 
centration dependence of D p , with RMS fractional errors 
in the range 4-18% (cf. Table HTJ) . 

The entirety of Lodge and collaborators' dataffi*^ on 
12-armed and 3-armed stars were also fit (separately for 
each arm number) to equation 1171 the joint stretched 
exponential in c, P, and M. Fitting parameters appear in 
IIIII Figure 23 shows the fit of D p for the 3-armed stars to 
the joint stretched exponential. As in the few other cases 
in which there are substantial differences between the 
fitted curve and measurements, the stretched-exponential 
form overestimates the concentration dependence of D p . 

Figure shows the outcomes of the joint fit to D p 
of the 12-armed stars. Except for the smallest (55kDa) 
star in the 1.3MDa matrix polymer, ea. 1 171 describes very 
well the entire dependence of D p on all three variables. 
The fit for / = 12 is markedly more outstanding than is 
the fit to the / = 3 stars, particularly with the 379 kDa 
pVME as the matrix polymer. The displayed curves rep- 
resent fits to all data points except for the 55kDa probe 
in the 1.3 MDa matrix polymer for c > 10 g/L. Including 
these 55kDa probe points in the fit raises the fractional 
RMS error from 14.6 to 19.5%, changes the three scaling 
exponents by 0.01 each, and elsewise has almost no effect 
on the fitted curves. 

Figures IT7I and E3%1 show that ea. 1171 and a single set of 
fitting parameters account well for the dependence of D p 
on c, P, and M for star polymers of given arm number in 
linear matrices. All three independent variables ranged 
over extensive domains: more than two orders of magni- 
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FIG. 47: D p of (a) 379 and (b) 1190 kDa 3-arm star 
polystyrenes in (from top to bottom, 140, 630, and 1300 kDa) 
polyvinylmethylethenortho-fluorotoluene as jointly fit to eq. 

EH 



tude in c, a factor of 30 in P, and an order of magnitude 
in M . The only failure in the fit occurred for the small- 
est probe in the largest matrix polymer (P/M rs 25) at 
elevated matrix concentrations. 

Martin2i*211 examined polystyrenes diffusing through 
polyvinylmethylether (as the matrix polymer) in toluene, 
with which the matrix is isorefractive. Polystyrenes had 
molecular weights of 50, 100, 420, and 900 kDa with 
M w /M n < 1.1. The polyvinylmethylether had from in- 
trinsic viscosity measurements a molecular weight ca. 110 
kDa and a 'fairly polydisperse' molecular weight distri- 
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FIG. 48: Dp of 12-arm star polystyrenes (from top to bottom, 
M w of 55, 467, 1110, and 1690 kDa) in (a) 140 kDa and (b) 
1.3 MDa polyvinylmethylether:ortho-fluorotoluene— , and the 
fit to the joint stretched exponential in c, P, and M. 



bution. Diffusion coefficients of the probe polymers were 
obtained with QELSS; Martin£i*2ii also determined the 
viscosities of the polymer solutions. 

Figure shows D p for each of the four probe poly- 
mers, as functions of matrix concentration. In Fig. 
149b . Dp for each probe polymer was separately fit to a 
stretched exponential in c, yielding the parameters given 
in Table [H] and the four solid lines seen in the Figure. 
Stretched exponential forms do an excellent job of de- 
scribing D p (c). In Fig. 149b. all data on the four probes 
was fit simultaneously to a single stretched exponential 
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FIG. 49: Probe diffusion of (top to bottom) 50, 100, 420, and 
900 kDa polystyrenes through 110 kDa polyvinylmethyl ether 
in toluene, based on data in Martin—*—. Solid lines are (a) 
fits of data for each probe and (b) fits of data for all probes 
to stretched exponential forms. 



in c and probe molecular weight P, obtaining fitting pa- 
rameters given in Table ITTT1 Martin only reports results 
for one M, making it impossible to evaluate the M- 
dependence of D p form his results. Plots of this stretched 
exponential, as functions of c at fixed P, give the solid 
lines of Fig. 149b . Agreement between the fitting function 
and experiment is good in the second Figure, though less 
good than in the first. 

Martin--' also measured the viscosity of his matrix 
polymer solutions. At lower polymer concentrations, es- 



FIG. 50: Dp of 342 kDa polymethylmethacrylate in 
polystyrene: thiophenol (data of Nemoto, et all— ) against 
polystyrene c concentration for polystyrene M of [top to 
bottom] 44, 186, 775, and 8420 kDa, and fits (a) for each 
M to a separate stretched exponential in c, and (b) to 
D P~ a exp(— ac v P 5 A'P) with D , a, and 7 as fitting param- 
eters. 



pecially for the lower-molecular- weight probe polymers, 
Martin 71 found that D p rj is nearly constant. At elevated 
concentrations, especially for the larger probe chains, 
Stokes-Einstein behavior ceases to obtain: D v r\ increases 
with increasing c, so that at large c and P the polymer 
chains diffuse faster than might have been expected from 
the macroscopic solution viscosity. 

In addition to the self-diffusion studies noted in the 
previous Section, Nemoto, et ali2£ used ultracentrifuga- 
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tion and quasielastic light scattering to measure sedi- 
mentation and probe diffusion coefficients of dilute poly- 
methylmethacrylates (as probe polymers) in isorefractive 
polystyrene: thiophcnol solutions. Polystyrene molec- 
ular weights were 43.9, 186, 775, and 8420 kDa; the 
PMMA had a molecular weight of 342 kDa. The reported 
M w /M n were in the range 1.10-1.17. 

Figure 150b shows Nemoto, et al.'s- measurements 
of D p of PMMA in solutions of each of the four 
polystyrenes. For each molecular weight of the matrix, 
D p (c) is described to good accuracy by a stretched ex- 
ponential in matrix concentration. Fitting parameters 
are given in Table ITT1 Figure ISTJb shows the same data 
with all measurements simultaneously fit to a stretched 
exponential in c and M, yielding parameters in Table 
IIIII Over nearly 200-fold variations in these variables, 
the single stretched exponential in c and M describes 
reasonably well the behavior of D p , with a 20% RMS 
fractional error. 

Nemoto, et ali^ report that D p /D a depends more 
strongly on c and P than does s/s . At large c, espe- 
cially at large M, D p /D a was found to be significantly 
less than s/s . Nemoto, et al. concluded that at elevated 
matrix concentrations and polymer molecular weights the 
sedimentation and self-diffusion behaviors of PMMA in 
polystyrene solutions are quite distinct. 

For two samples with the same D p /D a but very differ- 
ent matrix molecular weights (44, 8420 kDa), Nemoto, 
et al. also measured the shear viscosity r], finding that r\ 
differed 'by more than two orders of magnitude' between 
the two samples. Nemoto, et al. thus showed that D p is 
not governed by the shear viscosity of the matrix solu- 
tion. (The original paper did not specify which solution 
was the more viscous. Note that the comparison is being 
made at fixed D p /D Ql not at fixed c, so the correspon- 
dence is not self-evident.) 

Nemoto and collaborators^^ also used forced 
Rayleigh scattering to study the diffusion of probe 
polystyrenes through polystyrene:dibutylphthalate solu- 
tions. A first studyS focused on self-diffusion and tracer 
diffusion of labeled polystyrene through 40 wt% solu- 
tions of very long chains (M/P > 5) and very short 
chains (M/P < 0.2). Thirteen polystyrenes having 
2.8 < M w < 8420 kDa and M w /M n < 1.07 (except for 
chains larger than lMDa, for which M w /M n was in the 
range 1.09-1.17) were used in the studies. 

Nemoto, et al.'s dataS appears as Fig. [53 For D s 
and for D tr of short probe chains in solutions of long 
matrix molecules, a best-fit to the molecular weight de- 
pendence of D gives parameters seen in Table IIIII the 
P-dependence of the prefactor being forced rather than 
obtained from the fit. Because all data is at the same 
concentration, a concentration dependence was not ob- 
tained. Similarly, if we fit the data of Nemoto, et aliS 
on D p of long probe chains in short matrix chains to 
D Q P~ a exp(— aP 1 ), we find that D p gains its molecu- 
lar weight dependence almost entirely through the factor 
P~ a . The data are fit well with a = 0.52, in which case 
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FIG. 51: D s (filled points) and D p (open points) from Nemoto, 
et alj— , Tables 2 and 3. Dashed line marks systems with 
M <C P; solid lines represent a joint fit to D s and to systems 
with M/P > 5. 



the best-fit gives 7 « 0.03. 

In a separate paper, Nemoto, et ali^i used forced 
Rayleigh scattering to measure D p of probe polystyrenes 
in dibutylphthalate solutions of high-molecular-weight 
matrix polymers at 13 and 18 % matrix concentration. 
The probe polystyrenes were in the molecular weight 
range 6.1 < M w < 2890 kDa, with polydispersities 
M w /M n < 1.17, and generally < 1.09. Table IV shows 
the fits to this data. RMS fractional errors were 12-13% 
for single concentrations, but ca. 30% for the fit to both 
concentrations. It should be stressed that D p covers more 
than five orders of magnitude, so the errors are not large 
relative to the total range of D p . We find a significant 
dependence of D p on matrix as well as probe molecular 
weight. With only a few points for any particular P or 
M, it is difficult to present the data as a simple figure. 

Numasawa, et ali^ used QELSS to study the diffu- 
sion of tracer polystyrene chains through index- matched 
solutions of matrix polymethylmethacrylates in benzene. 
Polystyrene molecular weights were in the range 185-8420 
kDa, with M w /M n in the range 1.04-1.17. Polymethyl- 
methacrylates had molecular weights 850-4050 kDa, with 
M w /M n < 1.08, except for the 850 kDa polymer, for 
which M w /M n was 1.35. PMMA matrix concentrations 
were a half-dozen values in the range 0-36 g/L. The tracer 
polystyrenes were dilute in all solutions. In addition to 
measuring the probe diffusion coefficients, Numasawa, et 
ali^ also report the zero-shear viscosity of the matrix 
polymers and (on the basis of static light scattering) de- 
terminations of the radii of gyration of the probe poly- 
mers as a function of matrix concentration. 
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Numasawa, et al's results 85 appears as Fig. [""21 Fig- 
ure E3i shows D p of the 420 kDa and the 8.42 MDa 
polystyrene in solutions of two polymethylmethacrylates, 
as a function of matrix concentration. Figure l5"""b shows 
D p of five polystyrenes in each of four polymethyl- 
methacrylates, all at a matrix concentration near 37 g/L. 
One observes that D p decreases monotonically with in- 
creasing probe molecular weight and with increasing ma- 
trix concentration and molecular weight. 

All eight solid lines in both figures represent a simul- 
taneous fit of a stretched exponential in c, P, and M to 
all of the data in both Figures and to one additional data 
point. Fit parameters appear in Table ITTT1 It is impor- 
tant to emphasize that the same fitting parameters were 
used to generate all eight curves in both Figures. From 
Fig. 152b . the stretched exponential captures well the o 
dependence of D p , and the variation of that dependence 
with P and M. From Figure"5"2~b. at fixed c, the stretched 
exponential captures reasonably well the dependence of 
Dp on M at fixed P, but at fixed M and c does less well 
at capturing the dependence of D p on P. 

Nyden, et aliSi used PFGNMR to determine probe dif- 
fusion of monodisperse (M w /M n < 1.1) polyethylene ox- 
ides (molecular weights 10-963 kDa) diffusing through 
aqueous solutions of lOOkDa ethylhydroxyethylcellulosc. 
The authors studied 1% and 6% solutions and a 1% chem- 
ically cross-linked gel. Figure 15*31 shows D p /D Q of the 
polyethylene oxides in 1% of the matrix polymer as a 
function of their hydrodynamic radii. The solid line in 
the Figure represents a stretched exponential 

D p /Dpo = Di exp(-aR fi ) (19) 

in probe hydrodynamic radius R, where a and 5 are a 
scaling prefactor and exponent, D p q is the diffusion coef- 
ficient of the probe polymer in pure solvent, and D\ is the 
probe diffusion coefficient for a nominal R = polymer 
chain. The best-fit parameters were D\ = 0.33, a = 0.33, 
and f3 = 0.57; the RMS fractional error in the fit was 7%. 
Nyden, et al. 86 note R = KP a with a = 0.53 for probes 
having molecular weight P, so equation ["2"] is equivalent 
to a stretched exponential in P 7 with 7 « 0.30. This 
value of 7 is consistent with values for 7 found for other 
systems, as seen in Table ITTT1 

Nyden, et alm& also examined probe diffusion in 6% 
solutions of their matrix polymer. dD p /dM does not 
change monotonically with increasing M. However, at 
the larger M at which the anomalous behavior occurs, 
the PFGNMR echo decays are no longer simple exponen- 
tials, implying that diffusive behavior has become more 
complex. 

Pinder 8 ? reported from PFGNMR measurements the 
tracer diffusion coefficient of styrene and low-molecular- 
weight probe polystyrenes (P < 2470 Da) through solu- 
tions of deuterated polystyrenes (10.7 < M < 430 kDa) 
in theta (cyclohexane) and non-theta (CCI4) solvents. 
Dp was measured as a function of matrix polymer con- 
centration for matrix concentrations up to 300 g/L. The 
matrix polymer polydispersities were M w /M n < 1.14, 
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FIG. 52: Dp of polystyrenes in polymethylmethacrylate: ben- 
zene (a) as a function of matrix concentration, with P:M of 
[top to bottom] 420 kDa:4.05 MDa, 8.42 MDa: 1.95 MDa, 
and 8.42 MDa: 4.05 MDa, and (b) as a function of molecu- 
lar weight, for [top to bottom] 0.42, 1.26, 3.84, 5.48, and 8.42 
MDa polystyrene probes in 36.7 g/L polymethylmethacrylate. 
The lines all show the same best-fit stretched exponential in 
c, P, and M. 

with in most cases M w /M n < 1.08. 

Figure 154b shows the tracer diffusion coefficient of 
styrene in polystyrene solutions. The matrix molecular 
weights were 68 and 200 kDa in the non-theta solvent and 
68, 87, 200, and 430 kDa in the theta solvent. Within 
the scatter in the data — there are not a large number of 
data points for any particular matrix M — D p does not 
appear to depend on M. D p {c) was fit both to pure and 
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FIG. 53: Dp/Dpo from Nyden, et all— for polyethylene oxides 
with a range of molecular weights diffusing in 1% solutions of 
lOOkDa ethylhydroxyethylcellulose, plotted as a function of 
the hydrodynamic radii of the probes in pure water. 



stretched exponentials in c. As seen in Table [HI there 
appears to be very little improvement in the fit on al- 
lowing v 7^ 1, so the Figure shows the fits to the pure 
exponentials in c. 

Figure 154b shows the tracer diffusion coefficient of 
styrene and three low-molecular-weight polystyrenes 
(P < 2470 Da) though a range of high-molecular-weight 
polystyrenes (10.7 < M < 430kDa). For each probe, ma- 
trix molecular weights were varied by factors of 3-5 with 
no apparent significant effect on D p . In all cases, M/P 
was in the range 15-20 or larger. D p (c) for each probe was 
fit to a separate stretched exponential, giving excellent 
results reported in Table II. The merged results for all 
four probes were also fit jointly to a stretched exponen- 
tial in c and P. The computed D p for styrene monomer 
underestimates D p (c) by a nearly constant multiplica- 
tive factor. Because styrene is little larger than a sol- 
vent molecule, a second fit including only the polymer 
probes was made. Parameters for both fits are in Table 
IV. The Figure shows the second fit, to the three poly- 
meric probes. The second fit has a modestly better RMS 
fractional error, and finds in the limit of zero matrix con- 
centration that D p ~ p-0.52^ if {.jjg s ty r ene monomer is 
included in the fits, in the same c — > limit D p ~ p-0-68 
is obtained. 

Smith, et alpS used fluorescence recovery after pattern 
photobleaching to measure the diffusion of labeled 33.6 
kDa polypropylene oxide (PPO) chains through solutions 
of unlabeled 32kDa PPO chains dissolved in a melt of 1 
kDa PPO chains. The probe and solvent were relatively 
monodisperse (M w /M n = 1.1) while the matrix polymer 



I I I L 

10 20 50 100 200 

c (g/L) 




10 30 100 300 

c (g/L) 

FIG. 54: Dp, using data of Pindes^, of (a) styrene monomer 
through solutions of high- molecular- weight polystyrenes in cy- 
clohexane (theta solvent, filled points) and CCU (non-theta 
solvent, open points), and pure-exponential fits, and (b) [top 
to bottom] styrene monomer and 580, 1200, and 2470 Da 
polystyrene polymers in CCU and a high-molecular weight 
polystyrene matrix polymer with (lines) fits of the polymeric 
probes data to D p — D P~ a exp(— acP 1 ). 

was relatively polydisperse (M w /M n = 1.6). The matrix 
concentration was varied all the way from dilute solution 
up to the matrix melt. Smith, et al.'s data§£ are shown 
in Fig. EHJ The solid line in the Figure is a stretched- 
exponential fit using parameters in Table [H] A single 
stretched exponential with constant parameters describes 
the concentration dependence of D p all the way from 
dilute solution up to the matrix melt. 
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FIG. 55: D p of 33.6 kDa labeled polypropylene oxide chains 
through solutions of 32 kDa PPO chains in a 1 kDa PPO 
melt, using results of Smith, et alp—, Figure 1. The solid line 
represents a stretched-exponential fit. 



Tead and Kramei4S studied diffusion of 255 kDa 
deuterated polystyrene through solutions of large- 
molecular-weight (93, 255, and 20 000 kDa) polystyrenes 
dissolved in a low-molecular- weight (10 kDa) polystyrene 
melt. The matrix polymer volume fraction covered the 
full range < cf> < 1. The temperature was 150 C for 
the 93 and 255 kDa matrix polymers and 175 C for the 
20 000 kDa matrix polymer. Probe diffusion coefficients 
were obtained using forward recoil spectroscopy to mea- 
sure concentration profiles of the probe molecules as a 
function of time. As seen in Fig. [5BJ D p is represented 
well by a stretched exponential in the matrix polymer 
concentration c, no matter whether the matrix molecular 
weight is less than, equal to, or far larger than the probe 
molecular weight. The agreement with the stretched- 
exponential form is good from low concentrations of a 
short (93 kDa) chain out to high concentrations (melt) 
of a very large (20 000 kDa) chain. 

Tinland and Borsali 90 used fluorescence recovery af- 
ter photoblcaching and quasielastic light scattering to 
make independent measurements of the probe diffusion 
coefficient of 433 kDa dextran through solutions of 310 
kDa polyvinylpyrrolidone in water. PVP concentrations 
ranged from to 120 g/L. The polydispersity M w /M n is 
1.5 for the matrix polymer but ca. 1.9-1.95 for the probe 
chains. Except perhaps at the very highest concentra- 
tions studied, values of D p from the two techniques do 
not agree. In the following analysis, we use D p as ob- 
tained from FRAP, since the values from this technique 
do not rely so heavily on detailed model assumptions of 
the relationship between the scattering spectrum and the 



FIG. 56: D p of 255 kDa deuteropolystyrene in solutions of 
[top to bottom] 93, 250, and 20 000 kDa polystyrene in molten 
10 kDa polystyrene, at 150, 150, and 175 C, respectively, from 
Tead and Kramer 89 Figures 4 and 6. Solid lines are stretched- 
exponential fits using parameters given in Table II. For clarity, 
the 150 C data has been multiplied by a factor of 100. 

underlying diffusion coefficients. Tinland and Borsali's 
data appear in Fig. 1571 with data points from the orig- 
inal paper— and a smooth curve showing the stretched 
exponential in c generated from the fitting parameters in 
Table ITU Agreement between the measurements and the 
functional form is good at all polymer concentrations. 

Wheeler, ct al. 91 studied tracer diffusion of linear 
polystyrenes having molecular weights 65, 179, 422, and 
1050 kDa (with M w /M n < 1.1) through a 1.3 MDa 
polyvinylmethylether matrix polymer, M w /M n « 1.6, in 
orthofluorotoluene. D p was determined using quasielas- 
tic light scattering, which was possible because the 
polystyrenes were present at trace concentration while 
the matrix polymer and solvent are isorefractive. Matrix 
concentrations covered the range 1-100 g/L. For this ma- 
trix polymer, l/[rf\ — 2.2 g/L, so much but not all of the 
data is in the range c > c*. 

The data and corresponding stretched-exponential fits 
are in Figs. l55k and 158b . with fit parameters in Tables 
IITIand lTni respectively. Figure IS'Sk shows fits made sep- 
arately for each probe polymer. As in the other systems 
discussed above, the concentration dependence of D p for 
each probe is described extremely well by a stretched 
exponential in polymer concentration. Figure ISSb shows 
the outcome of a simultaneous fit of a single stretched ex- 
ponential in c and P to all data on all four probes. At the 
scale of the Figures, the curves in the two Figures for the 
179 and 422 kDa probes are coincident, while the curves 
for the 1050 kDa probes are very close. However, for the 
smallest probe and large matrix concentrations the joint 
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FIG. 57: D p of 433 kDa dextran diffusing through aqueous 
solutions of 310 kDa polyvinylpyrrolidone, based on data of 
Tinland and BorsaliSfi, Table 2 and Fi gure 3. Solid line is a 
stretched-exponential fit from parameters in Table ITTl 



stretched exponential substantially underestimates D p . 

Wheeler and Lodged used QELSS to measure 
the diffusion of linear polystyrenes through polyvinyl- 
methylether:orthofluorotoluene. The polystyrenes had 
molecular weights of 65, 179, 422 and 1050 kDa, 
with polydispersities M w /M n < 1.1. The polyvinyl- 
methylether samples had molecular weights of 140, 630, 
and 1300 kDa; concentrations of these matrix polymers 
ranged up to 300 g/L. Reference l92l represents a major 
extension of Ref. [9l| in the range of matrix concentra- 
tions, number of matrix molecular weights, and number 
of concentrations studied. For the three matrix polymers, 
c* was estimated at 11, 5.7, and 3.3 g/L, respectively 
(based on c* = 1.5/ [i]}), while c e was estimated at 50, 
12, and 6 g/L, respectively. The matrix polymers had 
M w /M n &l£. 

Wheeler and Lodge also used PFGNMR to measure 
D s q of orthofluorotoluene diffusing through the 1300 kDa 
polyvinylmethylether at matrix concentrations up to 300 
g/L. D s o fell by 62% over this concentration range. To 
good approximation, D s q(c) of the orthofluorotoluene is 
fit by a simple exponential D s o = 2.15 • 10~ 5 exp(— 2.98 ■ 
10 _3 c 1 ). Wheeler and Lodged and also Lodge, et alS 
used these data to modify D p of their polystyrene poly- 
mer probes to remove the concentration dependence of a 
nominal local friction £ = ksT / D s q. This local friction 
modification was here removed from the probe diffusion 
data of Refs. |9^ and [7^ before analyzing them further. 

Figures 1551 show Wheeler, et al's*£ data on their four 
probe polystyrenes in the 140 and 630 kDa polyvinyl- 
methylethers. Solid lines represent fits to stretched ex- 
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FIG. 58: Dp of polystyrenes (from top to bottom, M w of 
65, 179, 422, and 1050 kDa) diffusing through solutions of 
1.3 MDa polyvinylmethylether:ortho-fluorotoluene based on 
data of Wheeler, et all— , with solid lines showing fits to (a) 
separate stretched exponentials at each P, and (b) jointly to 
a stretched exponential in c and P. 



ponentials in c for individual probe:matrix pairs. The 
fits are uniformly extremely good, with RMS fractional 
errors in the range 2-12% and fitting parameters as seen 
in Table HD 

We also fit all of Wheeler, et al.'s dataSiiSS to a joint 
stretched exponential in c, P, and M, as seen in Figs. 1601 
The solid lines are the best-fit D p (c), plotted for fixed 
P and M, all lines being computed from a single set of 
parameters given in Table ILTTI In Fig. 160b . measurements 
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FIG. 59: D p of polystyrenes (from top to bottom, 
M w of 65, 179, 422, and 1050 kDa) diffusing through 
solutions of (a) 140 kDa and (b) 630 kDa polyvinyl- 
methylethenorthofluorotoluene based on data of Wheeler, et 
ali—, and (lines) separate stretched-exponential fits for each 
P : M combination. 

from Ref. |9^ were supplemented with the data of Ref. 0; 
the latter reference reported D p for all probes in 1300 
kDa polyvinylmcthylether solutions having c < 100 g/L. 

Over an order of magnitude in matrix molecular 
weight, a factor of 15 in probe molecular weight, and a 
factor of 300 in matrix concentration, the joint stretched 
exponential of eg 1171 represents reasonably well the joint 
dependence of D p on c, P, and M. The RMS fractional 
error in the fit is 25%. As noted above in the analysis 
for Ref. I^Js data, for the smallest probe in the 630 and 



1300 kDa matrix polymers (Fie. I60|) . the stretched ex- 
ponential form noticeably underpredicts D p . A similar 
issue arises for the 1050kDa probe in the 140 kDa ma- 
trix polvmer(fig. ItiOb) : at large c the predicted D p is too 
small, because the predicted curve does not bend quite 
sharply enough between small and large c. 



V. OTHER EXPERIMENTAL STUDIES 

In addition to the work reviewed above, a variety of 
other studies of polymer tracer diffusion appear in the 
literature. These are papers whose experimental foci are 
not the same as those of the papers examined in the pre- 
vious sections, and which are not amenable to analysis 
on the lines employed above. For example, a number 
of papers report scattering spectra of polymer: polymer: 
solvent mixtures in which neither polymeric species is di- 
lute, and therefore in which the spectral relaxation times 
do not correspond to the self-diffusion coefficient. 

Aven and Cohens^ measured the diffusion of dilute 
polystyrenes through 15 vol% solutions of polydimethyl- 
siloxane in tetrahydrofuran. The tracer diffusion coeffi- 
cients of the polystyrenes, the molecular weight depen- 
dence, and the initial linear slope of the concentration 
dependences were obtained with light scattering spec- 
troscopy. 

Borsali, et alMl studied QELSS spectra of 970 kDa 
polystyrene: 950 kDa polymethylmethacrylate: toluene, 
PMMA and toluene forming an isorefractive pair. Semi- 
quantitative experimental tests were made of the theo- 
retical work of Benmouna, et ali^i. When neither poly- 
mer was dilute, the observed spectrum was biexponen- 
tial. The mode amplitude ratio and relaxation times were 
within a factor of two of predictions of the Benmouna 
model—. The mutual diffusion coefficient of polystyrene 
in toluene and the cooperative diffusion coefficient of the 
mixture, both measured at the same total polymer con- 
centration, are equal to within 6%, also in agreement 
with the theory. 

Borsali, et alM^ extended their work on this ternary 
system with measurements at several large (weight frac- 
tion > 0.8) polystyrene concentrations and a range of 
total polymer concentrations, finding two relaxational 
modes in QELSS spectra. The diffusion coefficent associ- 
ated with the fast mode increased with increasing poly- 
mer concentration. The diffusion coefficient associated 
with the slow mode decreased markedly with increasing 
polymer concentration. Results were consistent with the 
Benmouna, et al^i model. 

Borsali, et al— i also studied mixtures of polystyrene 
and polydimethylsiloxane in tetrahydrofuran (which is 
isorefractive with polydimethylsiloxane) and in toluene, 
which is a zero average contrast solvent for these poly- 
mers and conditions. Measured spectral forms (one or 
two relaxations), relaxation times, and their concentra- 
tion dependences were consistent with the Benmouna 
model—. 
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In an extremely important series of experiments, 
Chang, et a\& report probe diffusion by polystyrene 
through polyvinylmcthylether: toluene, using both 
QELSS and FRS to determine diffusion coefficients. The 
two techniques differ markedly in the length scales they 
examine, scattering vectors q for the two techniques 
being I0 10 < q 2 < I0 11 cm' 2 for QELSS but only 
I0 7 < q 2 < I0 8 cm' 2 for FRS. The relaxation rates T 
were within experimental error proportional to q 2 over 
the full range I0 7 < q 2 < I0 11 cm -2 , confirming that 
QELSS does measure a simple translational diffusion co- 
efficient, even though the probe displacements sampled 
by QELSS are smaller than the diameter of the probe's 
correlation hole in the solution. Any claim that probe 
motion is more rapid, inside the probe's correlation hole, 
than it would be over larger distances, must therefore 
explain how Chang, et al.'s data are consistent with the 
claim. Chang, et ali2£ also demonstrate that the initial 
dependence of the measured diffusion coefficient on the 
probe concentration depends strongly on matrix concen- 
tration, the initial slope being substantially positive in 
the absence of matrix polymer and becoming significantly 
negative as polystyrene concentration is increased. 

B. Chu, et alii2i*i2£ report on the diffusion of 
polymethacrylate probes through a polystyrene ma- 
trix. The solvent was a mixture of toluene and a- 
chloronaphthalene. By adjusting the solvent composition 
and the temperature, it was possible to make a virtually 
exact match of the indices of refraction of the mixed sol- 
vent and the matrix, so that scattering arose only from 
the probe. By varying the solvent composition, the mo- 
tions of the matrix could separately be examined. A weak 
fast mode and a dominant slow mode were apparent, even 
at low angles, for probe chains much larger than the ma- 
trix chains. The slow mode had a weaker concentration 
dependence than does the solvent viscosity. 
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FIG. 60: Merged data from Figs. 0H3 and EH with (a) 140 
kDa, (b) 630 kDa, and (c) 1300 kDa matrix polymers, probe 
symbols as in the previous figures, and (solid line) fit of all 
data to a single stretched exponential in c, P, and M. Pa- 
rameters are in Table ITTT1 



S. Chu and collaborators** studied the motion of flu- 
orescently labelled DNA molecules in non-dilute solu- 
tion. They propose that they have confirmed the repta- 
tion model, based on their observations that (i) stretched 
DNA chains relax as though they were confined to tubes, 
(ii) D s scales as M~ 2 for nondilute probe chains, and (iii) 
a scaling law D s ~ c' 175 is consistent with their data. 

However, referring to these points seriatim: (i) Even 
within the reptation model a stretched chain relaxes from 
stretching primarily via its higher-order modes, revealing 
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only tube confinement on a short time scale, (ii) The 
scaling result D s ~ M~ 2 is a common feature of many 
models, is in no sense a signature of reptation^, and 
in any event is inconsistent with the moderniffi analysis 
of data on synthetic polymers, (iii) Ref. |98fs measure- 
ments of D s (c) are too scattered to demonstrate c~ 1,75 
behavior. Furthermore, if these data had demonstrated 
concentration scaling, they would be inconsistent with es- 
sentially all other data in the literature on the concentra- 
tion dependences of D s and D p of polymers in solution, 
as reviewed here, in that in almost all other systems D s 
and Dp have stretched-exponential, not power-law, con- 
centration dependences. 

Finally, as was not emphasized in the original paper, 
the data of Chu and collaborators*! 8 - conclusively reject 
the reptation model for their system. According to the 
reptation model, polymer chains only perform free diffu- 
sion after they escape from their tubes, which occurs at 
times longer than T<j. At times t < Td, polymer chains 
remain confined to their tubes and perform confined dif- 
fusion, so for t < Td a chain's mean-square dis plac ements 
satisfy (a; 2 ) ~ t x where x is 0.5 or less. Ref. |98| reports 
that their chains under their conditions remain within 
the hypothesized tube for a time w 1.2 or > 2 min- 
utes, depending on which of several methods is used to 
estimate Td- However, ref. l98l measured directly {x 2 (t)}. 
Figure 2 of ref. |9^ shows (a; 2 (_£}} ~ t 1 for times as short 
as Td/7 or less. That is, ref . |98| found that their chains 
were performing free, Brownian, non-constrained diffu- 
sion even even on time scales sufficiently short that under 
the reptation model the chains would have been confined 
to their tubes. This last finding is completely incom- 
patible with the reptation model, which requires that 
chains perform constrained curvilinear diffusion along 
their tubes so long as they remain inside their tubes. 
The data of Chu and collaborators^ 8 , thus indicates that 
a key aspect of the reptation picture, namely tube con- 
finement, is incorrect, at least in Chu, et al.'s system. 

Corrotto, et ali^fi^ performed static light scattering and 
QELSS measurements on bidisperse nondilutc mixtures 
of polystyrene in toluene, extracting fast and slow mode 
relaxations. 

Cosgrove and Griffiths^ used PFGNMR to study the 
diffusion of protonated polystyrenes through solutions of 
deuterated polystyrenes, varying the matrix concentra- 
tion, and the probe and matrix molecular weights. Data 
was obtained over limited ranges of D p on the depen- 
dence of Dp on each of these variables. D p generally de- 
clined with increasing matrix concentration and molecu- 
lar weight, and also declined with increasing probe molec- 
ular weight, except for very large probes in relatively di- 
lute solutions of small polymers, in which D p was nearly 
independent of matrix molecular weight. 

CottsiS^ report QELSS measurements on polystyrene: 
polyvinylmcthylether: toluene, proposing that radii of 
gyration, hydrodynamic radii, scaling exponents, and 
modes of diffusive motion can be measured systemati- 
cally for the dilute visible chains. This and other— early 



work on polystyrene: polyvinylmethylether: toluene 
flowered into the systematic studies of Lodge, Wheeler, 
and collaborators££'Z2'22i2ii22 on this system . 

Daivis and Pinde»i2& report QELSS studies of 
polystyrene: polyvinylmethylether mixtures dissolved in 
toluene or carbon tetrachloride. The first ternary mix- 
ture is nearly unique, in that the polymers are nearly 
compatible, toluene is a good solvent for both polymers, 
and toluene and PVME are isorefractive. The second 
ternary mixture differs from the first in that CCI4 is only 
a marginal solvent for polystyrene. In CCI4, spectra be- 
come nonexponential. The diffusion coefficient in nondi- 
lute CCI4 solutions is reduced by more than three-fold 
relative to diffusion by the same concentrations of the 
two polymers in toluene, but the probe radius of gyra- 
tion is only very slightly reduced by increasing the matrix 
concentration. 

Desbrieres and collaborators-^ applied QELSS to so- 
lutions of dextran and polyvinylpyrrolidone in water in 
solutions more concentrated than the overlap concentra- 
tion. Two modes whose properties are consistent with 
BcnmounaSi-type models of diffusion in ternary polymer 
solutions were observed. The mixtures have a phase sep- 
aration at elevated concentration. As polymer concentra- 
tions are increased towards the phase separation, a third 
slow mode whose various properties are consistent with 
the formation of aggregates was observed. 

Giebel and co- worker s 108 ! 109 studied QELSS spectra 
of polydimcthylsiloxane and polymethylmethacrylate in 
several solvents as a function of the relative concentra- 
tion of the two polymers. At the fixed total polymer 
concentration, the polymers were reasonably expected to 
be non-dilute. Comparison was made with theoretical 
results of Bcnmouna, et al—, with particular attention 
to the "zero average contrast" condition. A strong vari- 
ation of some spectral parameters with composition was 
described well by Benmouna—type models and a small 
number of free parameters. 

Jamil, et alii - report QELSS measurements on the 
diffusion of random-coil polystyrene probes through so- 
lutions of the rigid-rod polymer poly-(7-stearyl a,L- 
glutamate) in its isorefractive solvent toluene. The poly- 
mers are incompatible, poor miscibility of rods with ran- 
dom coils being identified by Jamil, et al—i as a signif- 
icant experimental challenge. Experimentally, the light 
scattering spectrum in these systems is due entirely to 
the polystyrene probes. 

Jamil, et al. 110 found that the dominant slow mode of 
the QELSS spectra is strongly concentration-dependent. 
At lower matrix concentration the diffusion coefficient 
corresponding to the slow mode increases with increas- 
ing probe concentration. At elevated matrix concen- 
tration this diffusion coefficient instead decreases with 
increasing probe concentration. Extrapolating this dif- 
fusion coefficient to zero probe concentration gives the 
tracer-diffusion coefficient of the probe polymer through 
the rod matrix. Jamil et al. reported fitting their six val- 
ues of D s to a stretched-exponential form. They found 
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a w 0.4 and v w 1.3. Jamil, et aliiiS interpret the find- 
ing v > 1 as arising from end-to-end aggregation of the 
matrix polymer; the matrix polymer increases its hydro- 
dynamic radius as its concentration is increased. 

Konak, et ali^ii report QELSS spectra of mixtures of 
polystyrene and polymethylmethacrylate in toluene. Nei- 
ther polymer was dilute. Comparison was made for a 
limited number of concentrations with theoretical mod- 
els arising from work of Benmouna, et al. 21 . Treating 
spectra as bimodal, the ratio of relaxation times was pre- 
dicted theoretically to better than 50%, but predictions 
of the mode amplitude ratio were often inexact by factors 
of 2 or 3. 

Marmonier and Lege»ii£ report extensive measure- 
ments using FRS on tracer diffusion of labelled 
polystyrenes through polystyrene in a good solvent. Un- 
fortunately, the reported data were modified by dividing 
them by an unreported concentration- dependent factor, 
namely the normalized concentration-dependent tracer 
diffusion coefficient of the free label in the same polymer 
solutions. It is therefore impossible to compare these 
measurements with other papers analyzed here. 

Numasawa, et al4i£ report scattering spectra of a se- 
ries of dilute polystyrenes in polymethylmethacrylate (in 
its isorefractive solvent benzene). The major focus was 
the scattering-vector dependence of the lincwidth T, es- 
pecially at larger scattering vectors q. At larger q, T/q 2 of 
the self-diffusive mode increases with increasing q. Mea- 
surements were only made at a single non-dilute matrix 
concentration, preventing further analysis of this data 
along the lines of this review. 

Russo, et aliii^ report self-diffusion coefficients for 
poly(7-benzyl-a,L-glutamate) in pyridine. The rodlike 
polymer has an isotropic-cholesteric liquid crystal phase 
transition with increasing concentration. Russo, et al. 
found that D s decreases with increasing polymer con- 
centration until the phase transition is reached. At the 
phase trnsition, D s increases adruptly; it then decreases 
again as the polymer concentration is further increased. 
D s (c) is qualitatively consistent with a stretched expo- 
nential form in the isotropic phase, but the number and 
spacing of points in the isotropic regime limits the accu- 
racy of a quantitative fit to this very interesting data. 

Scalettar, et alii^ used FRAP and FCS to study dif- 
fusion of phage A DNA solutions. By comparing sys- 
tems in which either few, or almost all, molecules were 
labelled, Scalettar, et al. were able to confirm the pre- 
diction of this author 13 that if very few macromolecules 
are labelled, fluorescence correlation spectroscopy de- 
termines their tracer diffusion coefficient, but if almost 
all macromolecules are labelled, fluorescence correlation 
spectroscopy determines their mutual diffusion coeffi- 
cient. 

Sun and Wan g 116 i 117 i 118 report a series of studies of 
polystyrene/polymethylmethacrylate mixtures (in ben- 
zene, dioxane, and toluene, respectively) using QELSS 
as the major experimental technique. Both polymers 
were in general nondilute. Neither polymer is isorefrac- 



tive with any of the solvents. The objective was to study 
the bimodal spectra that arise under these conditions and 
to show that the two relaxation times and the mode am- 
plitude ratio can be used to infer diffusion and cross- 
diffusion coefficients of the two components. Experi- 
mental series varied both the total polymer concentra- 
tion and the concentration ratio of the two components. 
The theoretical model predicts a biexponential spectrum; 
the experimental data was fit by a bimodal distribution 
of relaxation rates or by a sum of two Williams- Watts 
functions. The inferred self-diffusion coefficients of both 
species fall with increasing polymer concentration, but 
the concentration ranges that Sun and Wang studied are 
too narrow for further interpretation. 



VI. ANALYSIS 

The above sections have presented a detailed exami- 
nation of nearly the entirety of the published literature 
on polymer self-diffusion and probe diffusion in poly- 
mer solutions. The dependences of D s and D p on poly- 
mer concentration, probe molecular weight, and matrix 
molecular weight have been determined. Some features 
of this literature are incidental consequences of the chem- 
ical identity of the polymer being studied. The objective 
of this Section is to extract systematic behaviors from 
the above particular results. Here we ask: If we rise 
above particular features determined by the identity of 
the polymer under examination, what are the ideal fea- 
tures common to self- and probe-diffusion of all polymers 
in solution? 

In the following: First, the functional forms of the 
concentration and molecular weight dependences of the 
self- and probe diffusion coefficients are considered. Sec- 
ond, having found that D s and D p uniformly follow 
stretched exponentials in c, correlations of the stretchcd- 
exponential scaling parameters with other polymer prop- 
erties are examined. Third, for papers in which diffusion 
coefficients were reported for a series of homologous poly- 
mers, we examine a joint function of matrix concentra- 
tion and matrix and probe molecular weights and what it 
reveals about polymer diffusion. Fourth, we examine the 
few cases in which stretched-exponential behavior is not 
seen, or in which particular features of a given system 
clarify the systematic behavior of the phenomenological 
parameters used to describe D s and D p . Finally, other 
results implying a generalized phenomenology for aspects 
of diffusion behavior are examined. 

First, in the above I have reviewed virtually the en- 
tirety of the published literature on self-diffusion and 
probe diffusion of random-coil polymers in solution. As 
seen from the Figures in the preceding sections, the con- 
centration dependences of D s and D p are essentially al- 
ways described well by stretched exponentials (eq. Iltj|) in 
the matrix concentration c. 

Correspondingly, scaling (power-law) behavior is 
clearly rejected by almost the entire published litera- 
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ture on polymer self- and probe-diffusion. On a log-log 
plot, a stretched exponential appears as a smooth curve 
of monotonically varying slope. In contrast, on a log- 
log plot, scaling (power-law) behavior would appear as 
a straight line. Almost without exception, log-log plots 
of real measurements of D s (c) give smooth curves, not 
straight lines. Power laws could be fit to reported data, 
but in almost every case the power law would only pro- 
vide an accurate description of a tangent to the data over 
some narrow range of concentrations. 

An observation that experimental data is described 
very well by a particular mathematical form does not 
prove that the form in question is physically significant. 
In principle it may be the case that several different 
mathematical forms describe, to within the actual ex- 
perimental error, the same data. However, the observa- 
tion that experimental data is uniformly not described by 
some mathematical form, to well beyond experimental er- 
ror, is good evidence that models that predict the math- 
ematical form are inadequate. Power laws appear on log- 
log plots as straight lines. As seen from the above 70+ 
figures, straight lines are almost never observed in plots 
of D s or D p against c or other physical variables. Corre- 
spondingly, the concentration dependence of D s and D p 
almost never shows scaling (power-law) behavior. Scal- 
ing models that predict or assume for polymer solutions 
that D s and D p follow power laws in c, P, and/or M 
are very definitely rejected by almost the entirety of the 
published literature on D s and D p . 

Second, referring to eqs^l there are systematic cor- 
relations between scaling parameters a and v and the 
solution variables P and M. These correlations reflect 
the concentration and molecular weight dependences of 
D s and D p . In particular: 

The scaling pref actor a depends strongly on M . Fig- 
ure |^ shows the scaling prefactor a from measurements 
of the self-diffusion coefficient and fits to eq. 1161 as plot- 
ted against polymer molecular M. The Figure, based 
on Table I, shows almost all data on linear and star 
polymers, with concentrations in g/L. Results from the 
one systemSi with a large-concentration phase transition 
are omitted. M varies over nearly three orders of mag- 
nitude; a varies over almost four orders of magnitude. 
While there is substantial scatter, the figure is consistent 
a power- law correlation between a and M. The solid 
line in the Figure is a best fit to measurements on linear 
polymers. It shows a = 2.45 • 1CT 4 M 110 with M in kDa. 

It is also possible to examine the correlation of a with 
the size of the polymer, as reflected, e.g., by the diffusion 
coefficient D a . However, the variation in D a from system 
to system arises in part from differences in the measure- 
ment temperature and solvent viscosity. To eliminate 
these effects, D was used to compute a nominal chain 
hydrodynamic radius 



Here ks is Boltzmann's constant and T is the absolute 
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FIG. 61: Scaling pre-factor a for linear polymers (O) an d 
3-, 8-, and 18- armed star polymers (A,D,*) as listed Table 
I, plotted against polymer molecular weight M. Solid line is 
a power law a ~ M . 
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FIG. 62: Scaling pre-factor a for linear polymers, Table I, 
plotted against the bare hydrodynamic radius R, ea. 1201 Solid 
lines represent power laws with exponents of 3 and 4. 



temperature. The solvent viscosity 77 was taken from 
standard tables and interpolated as necessary to actual 
temperatures. In a few cases, reported experimental con- 
ditions do not permit an accurate conversion from D to 
R; these cases are not considered further. 

Figure plots a as a function of R. In Fig. [fJH the 
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FIG. 63: Scaling pre-factor a as a function of M for data 
from Refs. (O)^ for polymers in CCk, («) 46 , (U^M, 
(A)S£, (A)SiSSi with linear chains, with / = 3, (x) 62 

with / = 8, (H)fift with / = 18, (OA (♦)«, (VA for poly- 
mers in C 6 D 6 , (T)^, (®)2& (0)&, for linear polybuta- 
diene, (0)S for / = 3 polybutadiene, (>)32i for PEO in water, 
and (<)" for xanthan in water. Other details as in Fig. 1611 



solid lines are best fits for the data on linear chains (filled 
circles) toa~ R n for n of 3 or 4. If a ~ [77] as predicted 
by the Adler model 3 ^, then n = 3 is expected. The hy- 
drodynamic scaling model 3 ^^ instead predicts n = 4. 
The fit with n = 3 gives a modestly better fit to the 
data than does the n = 4 fit, but neither fit is visibly 
inconsistent with experiment. 

Figures IrTTl and Irj21 also present a of 3-, 8- and 18-armed 
stars, a of an 18-armed star polymer tends to be some- 
what smaller than a for a representative linear chain. 
However, a of star polymers almost always lies within the 
scatter in the values for a observed for the linear chains. 
a for three-armed stars includes results^ in which only a 
limited number of concentrations were studied for a given 
polymer. 

Figures and E21 emphasize polymer topology: one 
point style each for linear, 3-, 8-, and 18-armed chains. 
For some purposes, identifying the points by refer- 
ence (and, hence, by chemical system and experimental 
method) is more useful. Figures HkII and 1^41 give a against 
M and R, respectively, with points labelled by reference, 
this time including results of Tinland, et aliSi on a system 
with a large-c phase transition. The correlation between 
a and molecular weight is somewhat better that the cor- 
relation between a and the inferred - in most cases, not 
directly measured - chain radius. 

The scaling exponent v depends on M at low M . Fig- 
ure 1651 shows the dependence of v on polymer molecu- 



FIG. 64: Scaling pre-factor a as a function of R. 
details as in Fig. 1631 
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lar weight. The filled circles refer to linear polymers as 
studied in refsi 41 i 42 i 43 i 46 i 47 i 52 i 53 i 56 i 62 i 64 i 65 i 66 , but excludes 
systemsSi in which there is a phase transition at elevated 
polymer concentration. As seen from the Figure, for 
polymers larger than 250 kDa v approaches very closely 
to 0.5. For smaller polymers, v is substantially scattered, 
but increases with decreasing M . 

Sixteen systems in Table I, including linear and three- 
armed star polymers, had v forced to 1.00 during the 
fitting process. In these fits v = 1 was forced because 
the data would not support more free parameters, or be- 
cause the fit with v as a free parameter did not have a 
significantly better root mean square fractional error in 
the fit than did the fit with v = 1 forced. All but four of 
these fits refer to polymers with M < 200 kDa. The cir- 
cumstance that successful fits with v forced to unity are 
largely found with polymers of lower molecular weight is 
consistent with the interpretation that v increases toward 
1.0 at small M. These fits are represented in Figure IrTSI 
by the open circles. 

Figure IrjSI also shows the molecular weight dependence 
of v for star polymers in solution. The data on stars 
are substantially confined to smaller polymers (M < 250 
kDa). However, a trend in v with increasing M is appar- 
ent. At very small M, v is substantially scattered around 
1.0. For M > 200 kDa, v for star polymers appears from 
the few points to be trending toward 0.5. 

Third, several papera 42 i 67 i 71 i 73 i 74 i 75 i 77 i 78 i 79 i 80 i 82 i 83 i 84 i 85 i 87 i 91 i 92 
report D s and D p for a series of homologous polymers 
with different molecular weights. Each of these papers 
describes a study made by a given method using consis- 
tent operating conditions and data analysis procedures, 
thereby avoiding data scatter arising from any practical 
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FIG. 65: Scaling exponent v for polymers of Table 1 plotted 
against the polymer molecular weight M. A large-molecular- 
weight asymptote v m 0.5, and a low-molecular- weight in- 
crease in v toward v = 1 are apparent for linear polymers 
(filled circles) and for 3- (triangles), 8- (squares), and 18- 
(asterisks) armed star polymers. Open circles refer to linear 
polymers in which an adequate fit was obtained with v = 1 
forced during the fitting process. 

experimental issues. To use these studies to obtain 
information about the molecular weight dependences of 
D s and D p , data from each paper were separately fit to 
eq 1171 which is a joint stretched exponential in c, M , 
and P. (For D s , P and M are the same variable.) As 
seen Sections III and IV, eqElwith parameters in Table 
IIIII fits most of each data set well. 

Some studies were limited in the range of matrix con- 
centrations that was explored, either by not going to 
very high concentration or by not making measurements 
at low concentration. Theseii*£i22iSiS studies are ex- 
cluded here from further consideration. Some studies 
are subsets of other studies, for example cases in which 
measurements were first made at fixed M for a range of 
P, and then extended to cover a range of M. To avoid 
double-counting experiments, data sets that are subsets 
of other data are not considered further here. 

Section IV described fits of D p to eq El Much of 
the probe diffusion data is represented well by this joint 
stretched exponential in c, P, and M. However, eq 1171 
is less accurate when M and P differ substantially. In 
particular, for M/P > 3 or for P/M > 3, eqEl tends to 
overstate the decrease in D p with increasing c. Specific 
results of Section IV pertaining to the range of validity 
of earm include: 

Brown, et al&'s data for D p /D p q < 10 -3 show an ex- 
perimental D p that is smaller than the calculated value 
from a fit to eqEl Also, for P/M w 80, the measured 



Dp is larger than predicted. Hanley, et al.'s data™ and 
Kent, et al.'sZi data ( all follow well eqEI s P ana M 
dependences. Kim, et al.'s** extensive results on the ef- 
fects of varying P and M indicate that eq El generally 
works well for M/P < 3, but overstates the dependence 
of Dp on M for M/P > 3. From this data, it appears 
that as c is increased there is a narrowing of the range of 
values of M/P for which ea 1171 works well. 

For M/P never larger than 3 or so, Lodge and 
Wheeler^Z, found that the joint stretched exponential 
worked well for linear polystyrene chains and also worked 
well for / — 3 polystyrene stars, both diffusing through 
linear-chain matrices. For large / = 12 polystyrene 
stars in a short linear polyvinylmcthylcthcr, Lodge and 
Markland— 's data follow a joint stretched exponential, 
even for P/M sa 12. If the comparison is made by arm 
rather than total molecular weights, Lodge and Mark- 
land's data is confined to a region P a /M a < 2. Lodge, 
et al™ examined the diffusion of small / = 12 stars 
through solutions of a large linear polyvinylmcthylcthcr. 
For M/P > 3 and elevated concentrations, the measured 
Dp is larger than expected from the joint stretched expo- 
nential. Lodge, et al£2i£i also studied three-armed stars 
diffusing through linear chains; with increasing matrix 
molecular weight, the joint stretched exponential predicts 
too small a D p for M/P > 3. 

MartinSiiSi reported D p of linear polystyrenes diffus- 
ing through linear polyvinylmethylether. Equation 1171 
fits well (Fig. l49fl the joint data, but thoughout this data 
M/P < 2.2. Nemoto, et al™. studied diffusion of PMMA 
through polystyrene solutions; the joint fit (Fig. 1501) of 
all data to ea 1171 was more accurate for 2 > M/P > 0.5 
then it was outside this range. Numasawa, et al.'s data™ 
on polystyrene in polymethylmethacrylate:benzcnc, Fig. 
1521 is described reasonably well by eqEl However, for 
Numasawa, et al.'s data the joint stretched exponential 
clearly gives a better description of the concentration 
dependence of D p than it does of the M-dependence. 
Wheeler, et aliSiSS report on linear polystyrenes diffusing 
through polyvinylmethylethers. Fits to their measure- 
ments (Figs. I5§1 and 1601) show that ea 1171 underestimates 
Dp for M/P > 7 and for P/M > 3. 

The joint stretched exponential thus describes well a 
great number of measurements of D p in systems with 
0.3 < M/P < 3. Outside this range, eq El predicts a 
concentration dependence for D p that is stronger than 
that found experimentally. Corresponding to each set 
of measurements are a set of scaling exponents and pre- 
factors. To what extent do these parameters show uni- 
versal rather than system-specific behavior? In eq El 
data were parameterized as D a M~ a exp(— ac" P 1 M s ). 
For each system, the exponent v was treated as a sin- 
gle constant, the dependence of v on M at small M, 
discussed above, being suppressed. This suppression will 
reduce the quality of the fits when a very wide range of 
molecular weights is studied, as will be seen below in one 
system. 

From Table IIIII (i) a is almost always in the range 
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0.5 to 0.6. (ii) v is generally in the range 0.5-0.75. In 
two cases, v was found to be larger (0.86, 0.99) than this 
range; in one case it was smaller (0.43). In a few cases, 
v = 1 was forced during the fitting process, rather than 
being obtained from the fit. (iii) For linear chains, 7 is 
almost always in the range 0.25-0.3. There are outliers 
at 0.14, 0.19, and 0.43; each of these is associated with a 
system in which v was unusually small (0.43) or unusually 
large (0.86, 0.99). For star polymers, 7 is w 0.15. (iv) 
For linear polymers and three-armed stars, S is usually 
in the range 0.25-0.30. 

Fourth, in all of two cases, concentration and molec- 
ular weight dependences were not, or were not necessar- 
ily, described by stretched exponentials. First, Nemoto, 
et ali5£ studied the M-dependence of D s at large fixed 
c. Their data show a transition from a stretched- 
exponential to a power-law molecular weight dependence 
for M > 800 kDa. Second, Tao, et alM studied very 
large polymers in highly concentrated solutions. Their 
data are almost equally well-described by power laws and 
by stretched exponentials, as seen in Figs. 1131 and!14l For 
this data, a of eq. 1161 is nearly independent of M, con- 
trary to the behavior of a in other systems. The lack of 
an M-dependence in a, and the ability of a power-law 
to describe this data, are consistent with the proposal 60 
that polymer transport coefficients have different phe- 
nomenologies in the melt/near-melt domain and in the 
solution domain here. 

Fifth, there is evidence for other correlations. 

Several papers that report D s or D p also report the 
solution viscosity r\. In particular, MartmZi*£S gives not 
only Dp but also r\ of solutions of his 110 kDa polyvinyl- 
methylether matrix polymer. Martin reported for 50 kDa 
probe polystyrene that the product D v r\ increased dra- 
matically (up to 6-fold) with increasing c. For larger 
probe polymers (100 kDa, 420 kDa), the increase in D p r) 
was less dramatic; with 900 kDa probes, "...this product 
is very nearly independent of concentration." MartinSi 
interpreted the dramatic increase in D p r] with increasing 
matrix concentration as arising from a crossover from 
Stokes-Einsteinian to reptational diffusion, the crossover 
occurring more readily for P < M than for P ^> M. 

D p , t/ and probe radius of gyration R g for the same 
set of solutions were reported by Numasawa, et alS, 
who studied polystyrene diffusing through polymethyl- 
methacrylate:benzene. For a range of matrix concen- 
trations and matrix and probe molecular weights, Nu- 
masawa, et al. identified a regime in which D p rjR g is 
approximately constant. For small probe chains diffus- 
ing through larger matrix chains but not large probe 
chains diffusing through smaller matrix chains, the prod- 
uct D p 7]Rg increases markedly (up to 100-fold) with in- 
creasing c and M. Numasawa, et alp^ interpreted their 
data as showing Stokes-Einstein type diffusion and a 
crossover with increasing matrix molecular weight to a 
region in which reptational motion was dominant. 

Martinli and Numasawa, et al. 85 propose that the 
failure of the Stokes-Einstein equation, with D p r] in- 



creasing markedly with increasing c, is associated with 
a transition to rcptation dynamics for large polymer 
chains. However, the diffusion of rigid spherical probes 
through high-molccular-weight polymers shows a phe- 
nomenology highly similar to the phenomenology for Dr\ 
is observed^; namely Dr\ increases dramatically with in- 
creasing c at large M, even though there is no possibility 
that spheres have a transition permitting them to diffuse 
via reptation. 

Like effects do not prove that like causes are at work. 
The existence of non-rcptating spherical probe-polymer 
systems in which Drj increases markedly with increas- 
ing c does not prove that polymer systems showing the 
same phenomenology are not reptating. However, non- 
reptating sphere-probe systems can have a Drj that in- 
creases markedly at large matrix c. Therefore, the obser- 
vation that DpTj for probe polymers in a polymer solution 
increases markedly at large matrix c is not evidence that 
the polymer solution has had a transition to reptation 
dynamics. 

Xuexin, et ali^ 6 - report D p for linear and / = 18 
star polymers having very nearly the same D Q , i.e., very 
nearly the same hydrodynamic radius. For these two 
polymers, the scaling parameters a and v are also very 
nearly the same, consistent with but not proving the in- 
terpretation that a and v are determined by chain size 
and not by chain topology. 

Kent, et ali^i examined D p for a range of probe as 
well as matrix concentrations. For much of their data, 
Kent, et al.'s probes were concentrated enough that D p 
manifestly depended on the probe concentration c p as 
well as on the matrix concentration c. At low c, D p in- 
creased with increasing c p . At elevated c, D p decreased 
with increasing c p . These results provide a specific target 
for theoretical investigation. Just as good hydrodynamic 
models^ can predict quantitatively the concentration de- 
pendence of D of hard spheres, so also a good hydrody- 
namic model for polymers should be able to predict the 
dependence of dD p /dc p on matrix concentration. 



VII. CONCLUSIONS AND DISCUSSION 

In the above, virtually the entirety of the published 
literature on polymer self-diffusion and on the diffusion 
of chain probes in polymer solutions has been reviewed. 
Studies that determined how D s and D p depend on poly- 
mer concentration and molecular weight were systemat- 
ically re-analysed. Without exception the concentration 
dependences of D s and D p are described by stretched 
exponentials in polymer concentration. The measured 
molecular weight dependences also compare favorably in 
most cases with the elaborated stretched exponential, eq. 
IT7I Only when P S> M or M 3> P is there a deviation 
from ea. 1171 that deviation referring only to the molecular 
weight dependences. 

Contrarywise, almost without exception the experi- 
mental data is inconsistent with descriptions of concen- 
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tration or molecular weight dependences of the single- 
particle diffusion coefficient in terms of the power laws 
predicted or assumed by scaling models. The published 
data is sufficient to reject scaling model descriptions of 
polymer diffusion. 

The concentration dependence exp(— ac") of 
exponential-type models for D s is found to be valid for 
all ratios of the probe and matrix molecular weights. 
This validity of the stretched-exponential concentration 
dependence over a wide range of P/M is consistent with 
the mathematical structure of the renormalization-group 
derivation^ of these forms. Namely, the derivation 
begins with a low-concentration pseudovirial expansion 

D s = D {l + k 2 c + k 3 c 2 ) (21) 

for the concentration dependence of D Sl c here be- 
ing the matrix concentration, and then uses the 
Altenberger-Dahler positive-function renormalization 
group^S method to extend the series from lower to larger 
concentration. Over a modest range of P/M around 
P ss M, coefficients k 2 and k 3 with a simple dependence 
on P and M are adequate, leading to eqEI F° r P ^ M 
or P <C M, the simple calculations of k 2 and k% of Ref. 
l3ril are inadequate, because they do not include polymer 
internal modes, so the values of the fej change. However, 
for P 3> M orPcM the renormalization group method 
still receives as input a polynomial having the form of eq 
1211 so it will still generate as output a stretched expo- 
nential in c for the concentration dependence of D p . 

It is sometimes argued that the success of the stretched 
exponential form in describing D s (c) arises from a partic- 
ular flexibility of the stretched exponential, so that the 
successes shown in the Figures are accidental. Claims 
that D Q exp(— ac v ) is 'unusually flexible', relative to 
other functional forms, are inconsistent with basic math- 
ematics: The stretched exponential describes the con- 
centration dependence with three free parameters. The 
function is not singular for real c and positive v. There- 
fore, the region of function space spanned by the set of all 
stretched exponentials can be no larger than the region of 
function space spanned by any other function with three 
free parameters. Correspondingly, the functional form 
used here for D s (c) is no more flexible than other three- 
parameter forms, so its success in describing D s (c) can- 
not be ascribed to an unusual flexibility of the stretched 
exponential form. 

On the other hand, there is no significant region of any 
set of measurements in which concentration scaling is ob- 
served, other than as tangents to smooth curves of mono- 
tonically decreasing slope. It can always be asserted that 
corrections to scaling, of whatever basis, mean that the 
predicted power-law slopes are only observed over narrow 
ranges of concentration, or are only true asymptotically 
in the limit of large molecular weight: 

(A) Until the scaling models are refined sufficiently so 
as to predict the concentrations over which the hypothe- 
sized power-law slope should be approximately observed, 
it is impossible to tell whether the experimental data are 



consistent with a scaling prediction. It is entirely inade- 
quate to show that the experimental data has, in some re- 
gion, the slope predicted by a particular power law model, 
because the observed D s (c) someplace agrees with a pre- 
dicted c~ x behavior for every single positive x over a 
very wide range. If one is allowed to draw the asymp- 
tote where one chooses, the predicted x thus ceases to 
be falsifiable. On the other hand, a refined model might 
very well predict the slope for a range of concentrations 
in which the data did not have the predicted slope. (B) 
There is no indication in the experimental data that ex- 
ponential models cease to be adequate adequate as the 
polymer molecular weight is increased. 

On the other hand, as shown by the extremely thor- 
ough experiments of Tao, et alS, as one moves from 
polymer solutions toward the melt, one encounters a re- 
gion of large c in which scaling models are at least ap- 
proximately correct. Scaling models therefore may very 
well-this review does not examine this question in detail, 
and there is very little data, even though ilpS is very good 
data-be asymptotically valid in the near-melt regime. 
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TABLE I: Fits of the concentration dependence of the self-diffusion coefficient D s of polymers in solution to a stretched 
exponential D a exp(— ac"). The Table gives the fitting parameters, the percent root-mean-square fractional error %RMS, the 
polymer:solvent system, and the reference. Square brackets "[ ■ ■ ■ ]" denote parameters that were fixed rather than floated during 
the least-mean-squares fitting process. Do is in cm 2 /s or in units of Do- Concentrations are g/L, except * c in volume fraction 
units. Polymers include dex-dextran, pB-polybutadiene, HpBD-hydrogenated polybutadiene, pDMS-polydimethylsiloxane, 
PEO-polyethylene oxide, pl-polyisoprene, pS-polystyrene, and xanthan. Solvents include THF-tetrahydrofuran. 
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2.9 
8.1 
20 
8.6 
1.5 
3.6 
6.7 
17 



M(kDa) 
2.3 
6.5 
16 
21 
26 
75 
76 
90 
161 
227 
281 
6.5 
8.7 
29 
87 
10 
37 
179 
498 
1050 
32 
46 
105 
130 
360 
61 
92 
193 
216 
344 
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6300 
302 
70.8 
251 
302 
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pS 
pS 
pS 
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/ = 
/ = 
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/ = 
/ = 
/ = 
/ = 
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= 2 P B:CC1 4 
= 3 pB:CCl 4 
= 3 pB:CCl 4 
= 3 pB:CCl 4 
= 3 pB:CCl 4 
= 2 pB:CCl 4 
= 3 pB:CCl 4 
= 2 pB:CCl 4 
= 3 pB:CCl 4 
= 3 pB:CCl 4 
= 3 pB:CCl 4 
= 3 pB:CCl 4 
= 3 pB:CCl 4 
= 3 pB:CCl 4 
= 3 pB:CCl 4 
THF 
THF 
THF 
THF 
THF 
THF 
THF 
THF 
THF 
THF 

= 18 pLCCU 
= 18 pI:CCl 4 
= 18 pI:CCl 4 
= 18 pLCCU 
= 18 pI:CCl 4 
= 18 pI:CCl 4 
= 18 pLCCU 
= 2 pI:CCl 4 
= 2 pLCCU 
= 2 pLCCU 
= 2 pI:CCl 4 



pS:tol 



Refs 

63 



63 
63 
63 
63 
63 
63 
63 
63 
63 
63 
63 



64 



65 



65 



66 



66 
66 
66 
66 
66 
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TABLE II: Concentration dependence of the probe diffusion coefficient D p for molecular weight P probes in solutions of 
molecular weight M matrix polymers. The fits are to stretched exponentials D D exp(— ac v ) in the polymer concentration c. 
The Table gives the fitting parameters, the percent root-mean-square fractional fit error %RMS, the system, and the reference. 
Square brackets "[ • ■ ■ ]" denote parameters that were fixed rather than floated during the non-linear least squares fits. Dq 
is in cm 2 /s or in units of Do; concentrations are in g/L. Abbreviations include dex-dextran, "f = n" for an n-armed star 
polymer, hyal-hyaluronic acid, mr-methyl red, oFT-orthofluorotoluene, pl-polyisoprene, PMMA-polymethylmethacrylate, 
PPO-polypropylene oxide, pS-polystyrene, pVME-polyvinylmethylether, PF-pulsed field gradient nuclear magnetic resonance, 
QE-quasielastic light scattering spectroscopy, thiop-thiophenol, and tol-toluene. Parentheses indicate parameters from fits of 
limited accuracy. 
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0.64 
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0.26 
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0.57 


3.1 


65 


1300 


P S:pVME:oFT 


ZL 


5.42 


• nr 7 




0.21 
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"Various matrix molecular weights with M 3> P 
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TABLE III: Concentration and molecular weight dependences of D s and D p for molecular weight P probes in solutions of 
molecular weight M matrix polymers (for D s , one has P — M) at concentration c. The fits are to stretched exponentials 
D P~ a exp(— q?c"P 7 M' 5 ), with the percent root-mean-square fractional fit error %R, the material, and the reference. Molecular 
weights are in kDa; concentrations except as noted are in g/L. Square brackets "[ ■ ■ ■ ]" denote parameters that were fixed 
rather than floated. Abbreviations as in other Tables. 
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[0] 
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0.19 


[0] 


12 


b 


1300 


b 


/ = 2 P S:pVME:oFT 


77 


3.0 ■ 10~ 4 


0.54 


1.07 


10~ 3 
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0.61 
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0.11 


12 


b 


b 


180 


P S:pS:DBP 


84 


8.41 ■ 10 -6 


0.64 
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0.29 
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b 


b 


P S:pS:DBP 


84 


1.65 ■ 10~ 4 
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d 
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b 
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[0] 


14 


c 
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b 
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0.39 
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0.22 


25 


c 


b 


b 
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a Various, with M/P > 10 

'Various, see text. 

C A11 four probes, see text. 

^Excluding styrene monomer, see text. 

e Not all data points, see text. 

fp = M, self diffusion. 
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TABLE IV: Molecular weight dependence of the self and probe diffusion coefficients D s and D p for molecular weight P probes 
in solutions of matrix polymers at a fixed concentration c. The fits are to stretched exponentials D exp(— aM 1 ) in matrix 
molecular weight M. The Table gives the best-fit parameters, the percent root-mean-square fractional fit error %RMS, the 
system, and the reference. Square brackets "[ • • • ]" denote parameters that were fixed rather than floated. Abbreviations as 
per previous Tables, and DBP-dibutylphthalate. 
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Refs. 


8.54 • 1(T 8 


0.64 


0.55 


17 


1800 


a 


pS:pS:tol 


21 


2492P"- ' 5 


1.37 ■ 1(T 2 


0.47 


4.0 


P 


130 


P S:pS:dbp 


& 


1771p -o.5 


1.07 ■ 10~ 2 


[0.5] 


6.9 


P 


180 


pS:pS:dbp 





"Various, see text. 



TABLE V: Concentration and molecular weight dependences of the probe radius of gyration R g for molecular weight P 
probes in solutions of molecular weight M matrix polymers as functions of matrix concentration c. The fits are to stretched 
exponentials R g o exp(— ac"), with the percent root-mean-square fractional fit error %RMS, the materials, and the reference. 
Materials include EB-ethyl benzoate, PMMA-polymethylmethacrylate, pS-polystyrene. 



# 9 o(A) 


a 


V 


%RMS 


P(kDa) 


M(kDa) 


System 


Refs. 


390 


4.11 ■ 10" 3 


0.99 


2.2 


930 


1300 


pS: pMMA: 


EB 


74 


397 


8.06 ■ 10" 3 


0.77 


1.1 


930 


70 


pS: pMMA: 


EB 




[395] 


4.49 • 10" 3 


0.64 


1.9 


930 


7 


pS: pMMA: 


EB 


2* 



